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Abstract 

We study edge dynamics in the presence of interlayer tunneling, parallel 
magnetic field, and various types of disorder for two infinite sequences of 
quantum Hall states in symmetric bilayers. These sequences begin with the 
110 and 331 Halperin states and include their fractional descendants at lower 
filling factors; the former is easily realized experimentally while the latter is a 
candidate for the experimentally observed quantum Hall state at a total filling 
factor of 1/2 in bilayers. We discuss the experimentally interesting observables 
that involve just one chiral edge of the sample and the correlation functions 
needed for computing them. We present several methods for obtaining exact 
results in the presence of interactions and disorder which rely on the chiral 
character of the system. Of particular interest are our results on the 331 state 
which suggest that a time-resolved measurement at the edge can be used to 
discriminate between the 331 and Pfaffian scenarios for the observed quantum 
Hall state at filling factor 1/2 in realistic double-layer systems. 
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I. INTRODUCTION 



The dynamics of the edge modes in quantum Hall systems has been a subject of great 
interest for some years Its appeal is multifold. The low energy excitations of the 

ideal quantum Hall states that give rise to the plateau in the Hall resistance exist only 
at the edges. There is a deep connection between the structure of the bulk ground states 
and the "universal content" of the edge dynamics which is captured mathematically in a 
relation between 2+1 dimensional Chern-Simons theories and 1+1 dimensional conformal 
field theories [||||. This connection in turn implies a non-trivial charge dynamics at the edge 
which is now supported by experiments ||[7| . Finally, this connection allows the logic to be 
turned around in deducing new quantum Hall states from an analysis of possible conformal 
field theories . 

In this paper we investigate the edge dynamics of two infinite sequences of quantum 
Hall states in statistically symmetric bilayer systems in which we supplement the universal 
content by the inclusion of interlayer electron tunneling, an additional magnetic field parallel 
to the layers and, most importantly, disorder. The chief interest of this problem is that 
interlayer tunneling is strongly affected by the non-trivial charge dynamics and thus serves 
as an "internal probe" of the latter. In a previous publication ||, henceforth I, we had 
studied the problem of a clean system that gives rise to a chiral version of the sine-Gordon 
theory which is exactly soluble for the two infinite families of states for which interlayer 
tunneling is not irrelevant: these are the mm'n Halperin states with m = m' = n + 1 and 
m = m' = n + 2. The second of these families was shown to exhibit a remarkable trifurcation 
of charged excitations on the edge with the appearance of two Majorana fermions with 
dynamically generated velocities. In this paper we consider the additional effect of disorder 
and a non-zero temperature on the dynamics as well as the possibilities of modifying the 
tunneling by means of an interlayer magnetic field or a gated transfer of charge between the 
layers; the latter two procedures are essentially equivalent as we will see below. 

While the additional complications, especially the disorder, do not allow a complete 
solution in the sense of finding distributions of correlation functions in interlayer fields and 
at finite temperatures, the chiral character of the dynamics still allows us to make substantial 
progress in ways that should be of considerable interest to readers with a background in the 
physics of interacting, disordered systems. Consequently, we have included some amount 
of technical detail in the paper. In order not to lose sight of the principal physical results, 
especially those on the 331 state and the double layer Pfaffian state that are experimentally 
testable, this introduction is followed by a summary of the "useful content" of the paper. 
Readers primarily interested in the bottom line may wish to stop their perusal at its end. 

Before proceeding to that summary, a brief discussion of the observables and the relevant 
correlation functions is in order. As we will make more precise below, we study systems that 
possess one edge mode per layer so their single layer analogs are the Laughlin states at 
v = 1/m with m odd. In those single edge presents three natural observables fLOf . 



The first is the ground state expectation value of the edge current which can, in principle, 



do interesting things if the flux through the bulk is varied [[ll]]. The second is the tunneling 
density of states, computed from the one-electron Green's function, and the last is the edge 
mode velocity measured in a time resolved experiment done at the edge, which enters the 
retarded density-density correlation function. In real life, the first is not experimentally 
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relevant, while the third is not terribly interesting when there is one mode that is unaffected 
by disorder or temperature. The first has been experimentally investigated |||7j to great 
effect if not theoretical satisfaction, see Ref. [|12j and references therein. 

One of our central contentions is that the collective mode structure is interesting in 
double layer systems, even in the minimal cases where there is only one edge mode per 
layer. This was already clear in the clean cases considered previously, as in the trifurcation 
alluded to above, and is also the case in the more involved and realistic cases studied 
here. Consequently, the computation of retarded density correlators central to time and 
layer resolved measurements at the edges will be a central concern. In addition we will 
also compute one electron Green's functions needed for tunneling measurements but they 
will turn out to be essentially insensitive to the perturbations that we consider. We will 
not compute edge currents, although our results on partition functions in I can be easily 
extended to do so in clean systems, and extensions to the cases studied here are also feasible. 

We should note that some of this work has technical connections to earlier work on single 
layer systems Jl3|;[14[] with multiple edge modes but the details are different and one of our 
sequences, inclusive of the 331 Halperin state relevant to experiments, has no analog in 
single layer systems. 

The outline of the paper is as follows. We begin with a review of the edge theory of clean 
bilayer systems in the presence of uniform interlayer tunneling at the edge (Section |11 A| ) . 
Next we present our results for the effects of an interlayer magnetic field (Section |II B| ), dis- 
order (Section [II Cj ), and a finite temperature (Section [II D| ). The experimental consequences 
are discussed in Section |II E| , and the details of the calculations are presented in Section |T] 
and the appendices. 



II. SUMMARY OF RESULTS 

In this section we present a summary of our results. Readers interested in the details of 
the calculations will find them in Section fTJ. We begin with a review of the edge theory of 
clean bilayer systems. Next we consider the inclusion of a parallel magnetic field, disorder, 
and a finite temperature. We conclude with an experimental proposal for an edge measure- 
ment that could determine the bulk state responsible for the v — 1/2 plateau observed in 
bilayer systems. 



A. Review of Clean Bilayer Edge Theory 

The system under study consists of two parallel layers of 2DEGs in a strong perpendicular 
magnetic field. The geometry is sketched in Fig. [IJ Specifically, we are interested in the 
edge excitations of the Halperin states described by the mmn wavefunction, 
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a</3 a,f3 

where Zi a is the complex coordinate of electron a in layer i. The integer m determines the 
correlations within the layers and the integer n specifies the interlayer correlations. These 
states are incompressible and thus the gapless excitations of the system are confined to the 
droplet edges, which have length L and are parameterized by the coordinate x. 

The edge theory contains two chiral Bose fields, a charged mode, C , and a neutral mode 
0. We denote the velocities of these modes by u CjI1 , respectively. Excitations of the charged 
mode correspond to charge being added to the edge from the bulk, whereas excitations of 
the neutral mode correspond to a transfer of charge between the edges of the two layers. We 
restrict our discussion to states for which both edge modes move in the same direction, the 
"maximally chiral" case, by requiring m> n. 

In I it was shown that in the presence of interlayer single electron tunneling at the 
edge, the Hamiltonian of the edge theory separates into a free chiral boson Hamiltonian 
for the charged mode and a chiral sine-Gordon Hamiltonian (xSG) for the neutral mode. 
The chiral sine-Gordon Hamiltonian depends on the tunneling strength, denoted A, and the 
parameter (3 = \Jl{m — n), which sets the period of the interaction term. Since the neutral 
mode Hamiltonian depends only on m — n, the set of all maximally-chiral bilayer states 
can be divided into sequences labeled by the value of this difference. In particular, we will 
concentrate on the 110 sequence, which contains all states with m — n = 1, and the 331 
sequence, composed of states with m — n = 2. The tunneling perturbation is relevant, in a 
renormalization group (RG) sense, for the 110 sequence and marginal for the 331 sequence. 

The 110 and 331 sequences were solved exactly in I. In particular the single-electron 
Green's function 

g..( t ,x) = -i(T :^(t,x)::^(0,0):>, (2) 
and the two-point function of the density-fluctuation operator 
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iVij{t,x) = (T Pi (t,x) Pj (0,0)) - ( Pi (t,x))( Pj (0,0)), 



(3) 



where and P i are the electron annihilation and charge density operators on the edge of 
layer i, respectively, were computed exactly at zero temperature and L — > oo. We reproduce 
these results here, adopting a self-evident matrix notation. For the 110 sequence we have 



G(t,x) 
-iV(t,x) 



II cos(\x / v-n) + ia x sin(A:r/t> n ) 



\2n{x - v c t + ie t )} m -^2n(x - v n t + 



2(2m- 1) [2tt(x - v c t + ie t )} 2 



+ 



%-a x ) 



cos(2Ax/f n ) 



[2tt(x - v n t + ie t )Y' 



(4) 
(5) 



where 1L is the 2x2 unit matrix and a x is the standard Pauli matrix. For the 331 sequence 
we find 



Q{t,x) = 



-iV(t,x) = 



1 



[2n(x - v c t + ie t )\ 



1 



'11 



1 

"^2 
■f ^ 



2ti{x — v\t + ie t/ 
1 



+ 



2n(x — vit + iet) 



+ 



'1 



(7 



4(m - 1) [27r(a; - v c t + ie t )} 2 4 27r(a; - v x t + ie t )2n(x - w 2 t + ie t ) 



(6) 
(7) 



where we have defined the velocities i>i,2 = f n =t \/ir, and e t = sgn(t)a, where a is a short 
distance cutoff. Each of these functions contains a part arising from the charged mode and a 
part from the neutral mode. For the single-electron Green's functions the contributions from 
each mode are combined multiplicatively, while for the density-density correlation function 
they are combined additively. 

In addition to these time-ordered correlation functions, in later sections we will be inter- 
ested in the corresponding retarded functions which govern physical response measurements 
at the edge. The density response function is 



V R (t,x) 



2tt 



(l+(T a 

2(2m- 



-5'(x — v c t) + 



cos(2Ax/f n )<5'(a; — v Q t) 



(8) 



for the 110 sequence, and 



V H (t,x) = 



2tt 



o (x - V c t) + 



a 



4(m - 1) 



4(i;i - v 2 )x 



{v 2 5(x - v 2 t) - vi5(x - vit)} 



(9) 



for the 331 sequence, where the prime denotes differentiation with respect to the argument. 

We see that for the 110 sequence relevant tunneling produces spatial oscillations in the 
correlation functions, while for the 331 sequence marginal tunneling leads to two velocities 
(1*1,2) in the neutral mode sector, and hence a total of three velocities for the system as a 
whole. Note that even at zero temperature and in the absence of disorder the signal from 
the neutral mode in the density-density response function decays with distance because of 
tunneling. In the following section we shall investigate how these correlation functions are 
modified by various perturbations. 
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B. Parallel Field 



We first discuss the effects of an interlayer magnetic field. If we take the z-axis along the 
normal to the layers and recall that the x-axis is along the edges, we consider an additional 
magnetic field along the ?/-axis: B = By. The edge Hamiltonian in the presence of a 
parallel field depends on the parameter Y = v n Bd(3/2, where d is the layer separation. We 
remark that the effect of the interplane magnetic field considered here is distinct from the 
simple decrease in A caused by the reduction of the interlayer tunneling matrix element. As 
noted by Chalker and Sondhi, these effects can be distinguished by studying large aspect- 
ratio samples [T5|]. The results here also apply to the case where we introduce an electric 



potential difference between the layers instead of an interplane magnetic field. 



1. 110 Sequence 

For the states in the 110 sequence we find that the spectrum of the edge theory in 
the presence of an interlayer magnetic field can be obtained from the spectrum with zero 
interlayer field via the replacement 



A i-^ A' = y A 2 + T 2 /2, (10) 

in particular we see that the number of velocities is unchanged. The two-point function of 
the density-fluctuation operator (H) is 

-W{t,x) = — 7T" 1 7 /. , 19 + 



2 (2m- 1) [2tt(x- v c t + ie t )} 2 2(A') 2 
(11 -a x ) fr 2 x2 (2X'x\\ 

x 7 7 — " TFTTS hA 2 COS >. (11) 

Note that in the absence of tunneling (A = 0), the correlation function is unaffected by the 
parallel magnetic field, i.e., it is independent of T, as expected. For non-zero tunneling, the 
addition of an interlayer magnetic field increases the frequency and decreases the amplitude 
of the spatial oscillations in the density-density correlation function. One can show that the 
effect of the parallel field on the single-electron Green's function is similar to its effect on 
the density two-point function, see Eq. 



2. 331 Sequence 

For the states in the 331 sequence we find that in the presence of an interlayer magnetic 
field the spectrum of the neutral mode portion of the edge theory is 

Hb = y ^s{k) :a\.ak-, where (12) 



e(k) = v n k + sgn(rW(Afc/vr) 2 + T 
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where the are canonical Fermi operators. Note that in the limit of vanishing parallel 
magnetic field, T — > + , the energy dispersion (for A > 0) becomes e{k) = (v n + \/ir)k = V\k 
for k > and e(k) = (v Q — X/ir)k = v 2 k for k < 0, which is the spectrum of two right-moving 
Majorana fermions with split velocities. (In the limit r — > 0~, or for A < 0, we get a similar 
dispersion with v\ and v 2 interchanged, but this does not alter the excitation spectrum of 
the Hamiltonian) . For any non-zero interplane field we find that the dispersion develops 
some curvature, which corresponds to the two Majorana species being mixed at distances 
large compared with A/7rr. The dispersion is sketched in Fig. [| 




k 



FIG. 2. The solid line is the dispersion e{k) plotted for T > 0. The dotted line is e = v 2 k and 
the dashed line is e = v\k. Note that e(k) asymptotes v\k as k — > oo and v 2 k as k — > — oo. 



The correlation functions for the 331 sequence with a parallel field can be reduced to 
quadrature. If we define the quantities r = sgn(t)(t> n X — t), 

X = k(u) = y/^i^r 2 + (A/tt) 2 ^ 2 , (13) 

and the function 

v ' ; Jo 2vr k(u) V ; 

then the single-electron Green's function and the density-density correlation function can 
be expressed in terms of P(r, X) and its derivatives. We find 
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G(t,x) 



-isgn(t) 



m— 1 



7T / 



wif 2 [27r(x - v c t + ie t )] 
+ sgn(t) [ -1 - ^a* ] -P r (r, X) - u^ 2 rV 2 P(r, X) 

V 7T / 7T 



-W(t,x) 



4(m - 1) [2tt(o; - w c t + ze 4 )] 2 



;a - ^ 



p2 



A 5 



7T 



p; 2 + ^ 2 r 2 p 2 



(15) 
(16) 



where the subscripts on P denote partial differentiation. 

In Fig. [I we plot the real and imaginary parts of the neutral mode part of Qw at fixed 
X = 10 as a function of — r for v n = 1, T = 0.5, and \/n = 0.1. The corresponding plot 
for the neutral mode part of Q12 is given in Fig. ^. The singularities visible in these plots 
occur at the points t = x/v\ and t = x/v^, corresponding to propagation at the Majorana 
velocities. We see that the parallel field does not destroy the velocity-split structure of 
the Green's function, which is somewhat remarkable. The parallel field is an RG relevant 
perturbation which modifies the low-energy, long-wavelength properties of the system, see 
Fig. 0. Nevertheless, the singularities at v\ and v 2 are completely unchanged, see Eq. (|166|) 
below, since they arise from integrations over all frequencies. In Appendix || we will discuss 
the opposite case: an RG irrelevant perturbation which does modify some features of the 
Green's function. At times corresponding to propagation velocities between v\ and f 2 we 
find oscillatory behavior not present at Y — 0. 




-2 - 



X/w=o.i,x=io, r=o.5 



Re(G) 11 
\m(G) u 



-2 



-T=t-v„X 

FIG. 3. The real (solid line) and imaginary (dashed line) parts of the neutral mode factor in 
Qn (74) plotted as a function of — r for v n = 1, X/tt = 0.1, X = 10, and T = 0.5. 
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FIG. 4. The real (solid line) and imaginary (dashed line) parts of the neutral mode factor in 
G12 ( [74] ) plotted as a function of — r for v n = 1, \/ir = 0.1, X = 10, and T = 0.5. 



C. Disorder 

We now consider the effects of disorder on the bilayer system. Our primary interest is 
how the novel features of the correlation functions in the presence of tunneling, i.e., spatial 
oscillations in the 110 sequence and the splitting of the neutral-mode Majorana velocities in 
the 331 sequence, are modified by disorder. Note that the quantities of physical interest are 
the correlation functions in a typical realization of disorder, whereas the readily calculable 
quantities are disorder-averaged correlation functions. The typical and average quantities 
may have very different behavior, and we will discuss such differences at various points. 

In addition to the possibility of disorder in the tunneling amplitude, A, we consider 
random scalar potentials, which couple to the edge charge densities in each layer. 

From a perturbative RG analysis we find that disorder in A is relevant for the 110 sequence 
and irrelevant for the 331 sequence, while the random scalar potentials, £i(x), are relevant 
for both sequences. 



1. 110 Sequence 



We consider a tunneling amplitude which has mean A and variance A\, and a disor- 
dered scalar potential with zero mean and variance Ag. We find that the disorder-averaged, 
retarded density response function is 



T> R (t,x,x') 



2tt 



1 + ^') r // / \ (1 -a 31 ) P „ 

o [x — x — v c t) H o [x — x — v n t) 



x e 



-2|x-x'|(A A +A 5 /4)/u 



2(2m - 1, 

f 2\x — x'lA 



cos 



+ 



A, 



4v n X 



^ sin 



2\x — x'\\ 



(17) 



9 



where 



A 



The disorder in the tunneling amplitude (Ax) produces an exponential decay with distance 
in the neutral mode part of the disorder-averaged density response function. The random 
scalar potential (Af) has a similar effect, and in addition it produces a shift in the frequency 
of the spatial oscillations ([18]). Using the fact that the neutral-mode part of P in a given 
sample can be expressed in terms of products of single-particle Green's functions, whose 
absolute squares are long-ranged (i.e., algebraically decaying) we can conclude that in a given 
sample the density response function has the structure of the disorder-averaged quantity 
([H]), without the exponential decay in space of the neutral mode piece. 



2. 331 Sequence 



Above we remarked that for the 331 sequence only disorder in the scalar potential terms 
was a non-irrelevant perturbation. We therefore consider only a disordered scalar potential 
with zero mean and variance A^. We find that the disorder-averaged, retarded density 
response function at T = is: 



T>*(t, x, x') = JL±^-5'(a; - x >- VB t) + ^—^D*(t, x, x') 

Z7T Aim — 1 4 



(19) 



where the neutral mode contribution is 



D R (t,x,0) 



-A^x/viV2 




— — {v\5(x - fit) - v 2 5(x - v 2 t)\ + - — —I 

(vi - v 2 )x 2v x v 2 



A, 



A/ 



7T 



71 T 



A/tt 



(20) 



Here z = (t — x/vi)(x/v 2 — t), I n are Bessel functions of imaginary argument, and to = 
(x/v\ + x/v 2 )/2 = v n x/viv 2 is the mean arrival time. 

Comparing this expression to the result for the neutral mode in a clean system, the second 
term in Eq. (gj), we find that the delta-function peaks at the velocities V\ and v 2 remain sharp 
in the presence of disorder, but their amplitudes acquire an additional exponential decay. 



In the disorder-averaged result, (|20|) , there is also a signal centered around the mean arrival 
time to. If we write r — to — t, then in the limit of large distances: x/v\v 2 ^1/ A^, for times 
near the mean arrival time, (tt/X)t x/v\v 2 , the central signal is asymptotically 



D'Ht.r.O) % \j^V ) cxi> 



A^v x v 2 ^ 



2x(A/tt) 



(21) 



This asymptotic form is similar to a result obtained by Wen for the case of two v = 1 
edges with unequal velocities 113]. The reason for this similarity is that both problems are 



formally equivalent to a spin in a random magnetic field which undergoes diffusion on the 
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SU(2) group manifold, see Appendix [A| The term in Eq. (pip decays algebraically with 
distance. Therefore, while the signal in D R at the extremal velocities (1*1,2) is exponentially 
suppressed by the disorder, there is an additional signal with velocity V\v 2 /v n which only 
falls off algebraically. 

To determine the behavior of D R in a given realization of disorder we have adopted 
several approaches. First, as in our analysis of the 110 sequence, we can use the fact that 
D R can be expressed in terms of single-particle Green's functions whose second moments 
we can evaluate. Second, we have found that D R (t,x,x') exhibits an exact antisymmetry 
about the point t = t in each realization of disorder: 

D R (t + t, x, x') = -D R (t - t, x, x'). (22) 

This is an interesting result because it is an exact dynamical symmetry (i.e., it is a relation 
between correlation functions at different times) in a system with an arbitrary potential 
£(x). Finally, we have calculated the behavior of the correlation functions for two simple 
potentials, the case of a uniform potential £(x) = const, and the case of isolated delta- 
function impurities £(x) = J2 m <lmfi(x — y m ). All of these results, discussed in detail in 
Section |111 Q lead us to the conclusion that D R (t,x,x') in a typical configuration is given 



by an expression similar to Eq. ( |20| ) for D R (t, x, x'), with the principal value factor replaced 
by a function f(t, x, x'), which is a rapidly fluctuating function of time, antisymmetric about 
the point t = to, and whose amplitude grows as t approaches to- These conclusions about 
the behavior of D R in a given sample have been verified by numerical simulations which are 
discussed in Section 111 F. 



D. Finite Temperature Effects 

We now briefly consider the effects of a finite temperature T = 1/(3. For a single chiral 
edge mode we know that at zero temperature V R (t,x) oc 6{t)8'(x — vt). A straightforward 
calculation shows that this form is actually temperature independent. 

Recall that for the 110 sequence the retarded density- density correlation function is a 
sum of terms of the form 6(t)5'(x — vt) multiplied by a function independent of t. This can 
be seen for the clean system in Eq. (^j), and for the disordered system in Eq. flT7|). We can 
therefore conclude that T> R for the 110 sequence is temperature independent even in the 
presence of a non-zero scalar potential. 

The situation is different for the 331 sequence. While the term in T> R from the charged 
mode is temperature independent for the reasons given above, the neutral mode term is not. 
At a finite temperature one finds for the neutral mode in a clean system: 

D R (t, x) = - JpL_— _ [v 2 5(x - v 2 t) - Vl S(x - Vl t)] . (23) 

2(3vxV 2 sinh [ti(v 1 — v 2 )x/ (3v\V 2 \ 

Comparing this with the zero temperature result, the second term in Eq. @, we find that 
the neutral mode term, which decays as 1/x at T = 0, decays exponentially at T > 0. 
In the presence of disorder one finds that the finite-temperature form of D R is 
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Z7T 



{ir/(3viv 2 ) 



sinh[7r(fi — v 2 )x/ (3v\v 2 ] 



{t>i<5(x - v\t) - v 2 5(x - v 2 t)} 



+ Az6(z)V 



{ir/Pv-LVz) 



' x/viv 2 



sinh[2vr(t -t)/(5] 



h 



A/tt 



71 T 



A/tt 



(24) 



This result was obtained using the formalism given in Appendix [A]. 

Comparing this with the zero temperature result in Eq. (|2~0l) we see that the amplitudes 
of the delta-functions at the extremal velocities acquire an additional exponential decay 
because of the finite temperature. However, the replacement 



V 



t Q -t 



V 



(27T//3) 



sinh[27r(t -*)//?]. 



(25) 



indicates that the structure in D R centered on the mean arrival time sharpens at a finite 
temperature. 



E. Experimental Ramifications 

One of the primary results of the previous sections is the unusual structure of the retarded 
density-density correlation function for the 331 sequence. The first state in the 331 sequence, 
the state 331 itself, has a filling factor of 1/4 per layer, for a total filling factor of 1/2. A 
plateau in the Hall conductance has been experimentally observed at v = 1/2 in bilayer 
systems jl6|| . Another candidate state which has been proposed to explain this plateau 



is the Pfaffian state fT7|. Standard experimental probes of the edge states, such as the 



non-linear I—V characteristic, cannot be used to distinguish the 331 from the Pfaffian state 



since both states give the same power law exponent |TB[. In this section, we argue that the 
retarded density-density correlation function of the Pfaffian is sufficiently different from that 
predicted for the 331 state that, even in the presence of a finite temperature and disorder, 
a measurement of this correlation function at the edge could distinguish between these two 
bulk states. 

For the Pfaffian edge theory we find that the retarded density response function at a 
finite temperature and in the presence of disorder is 

V&(t, x) = -^-8{t)6'{x - V)- (26) 

47T 

We see that there is only a single velocity present. 

Recall from Section [II A , Eq. @, that for the 331 edge there are three velocities present 



in the clean system at zero temperature, one for the charged mode (v c ) and two for the 
neutral mode (^1,2), whose splitting is due to tunneling. The signal at the two neutral 
mode velocities decays as 1/x at T = in the clean system. In Section |11 C] we saw that a 
disordered scalar potential suppresses the signal at the extremal velocities by a factor which 
decays exponentially with distance. However, we found that there is a broad signal in the 
neutral mode, centered on a velocity distinct from the charged mode velocity, which decays 
only algebraically. In Section |IID| we saw that a finite temperature actually sharpens the 
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signal centered on the mean arrival time. From the results and discussion in Sections [TIC 



and [II D| , we expect that for the 331 edge in a given realization of disorder, the structure of 
the retarded finite-temperature density- density correlation function is: 



f 31 {t,x,x') « 0(t)(U + a x )S'(x - v c t) 



f{t,x,x',P) 



£(x 



exp 



x 



(t - 1 ) 2 
£(x — x') 



(27) 



where the function f(t,x,x',/3) is the finite-temperature version of the function introduced 
at the end of Section [11 (J| , and £ is a length scale set by the potential. We expect f(t, x, x', /?) 
to have the same properties as the zero-temperature function, and to have its support more 
strongly concentrated near t = t Q as the temperature increases. 

The experiment we propose involves creating a density disturbance at one point along 
the edge of the bilayer system and measuring the signal some distance downstream. The 
experimental geometry is sketched in Fig. |5]. Basic linear response theory states that if the 
density disturbance is produced in layer j via an external potential V ex (t,x), and measured 
in layer i then the signal is 



(pi(t,x)) c 



dt'dx Vj^At — t' , x, x')V e x(t' , x). 



(28) 



If the external potential is turned on at a point, i.e., V ex (t',x') = 5(x')6(t'), then the mea- 
sured signal is 



(Pi(t,x)) m = f dt'V^(t',x,0). 

J — oo 



(29) 



AAA 




B 




FIG. 5. The experimental geometry showing the quantum Hall bilayer in a magnetic field B 
with two spatially separated contacts. A density disturbance is produced at one electrode (Mn) 
and measured at the other electrode (V OVL t)- 



For the Pfaffian state one would see a single sharp signal, see Eq. ([26|). In contrast, for the 
331 state, in addition to a sharp signal from the charged mode there would be a second signal 
from the neural mode. To illustrate the neutral mode signal one would expect in this case we 
have performed numerical simulations. A typical trace, computed at zero temperature for 
v n = 1, \/n = 0.1, X = 10, Ag = 0.1 with stepwise-constant disorder potential with N x = 2 7 
values is shown in Fig. [] with a solid line. As expected, it is a rapidly fluctuating function, 
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but it exhibits an exact symmetry about the mean arrival time. This symmetry follows from 
the antisymmetry of D R fl22|) and the time integration (|29"D . Although the signal is very 
noisy, a measurement with finite resolution (dashed line) produces a curve which does not 
average to zero. The amplitude of the smoothened signal is maximal at the mean arrival 
time as was predicted in Section |11 C| . In Fig. |7| we show the result of numerically averaging 
over 1600 impurity configurations (solid line), as well as the analytic average (dashed line) 
evaluated using Eq. (P0"|). The two curves are in excellent agreement. The finite width of 
the neutral mode signal is a novel feature of the 331 state; in the 110 sequence the neutral 
mode propagates with only one velocity. 

In Fig. [8] we show a similar trace, but with the disorder stronger by an order of magnitude, 
A^ = 1.0. Note that the amplitude of the signal near the extremal arrival times r = ±1 
is suppressed relative to the case with a smaller disorder strength, but the signal near the 
mean arrival time (r = 0) is not. In Fig. ^ we show the result for A^ = 0.1 at a high 
temperature T/A^ = 2000. The signal in a given realization of disorder (thin solid line) is 
as noisy as in the zero-temperature case (Fig. [(J), however the amplitude of the smoothened 
signal (dashed line) is down by roughly an order of magnitude. At this high temperature 
the disorder-averaged result (thick solid line) is essentially zero everywhere except very close 
to the mean arrival time. 




- 1 "0-5 o -t=m>, x 0-5 1 



FIG. 6. The neutral mode contribution to the integrated T>i 2 (^) at T = for v n = 1, 
A/-7T = 0.1, X = 10, and = 0.1. The horizontal axis is time measured from the mean arrival 
time. The solid line is for a given realization of disorder, and the dashed line assumes a measurement 
with a finite time resolution. 
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-0.4 



A/jc=0.1,X=10, Ap=0.1 
_j j 



average (yv=1600) 
analytic average 



-0.5 



0.5 



FIG. 7. The neutral mode contribution to the integrated P(| (P^l) at T = for -u n = 1, 
A/-7T = 0.1, X = 10, and A^ = 0.1. The solid line is the numerical average over 1600 impurity 
configurations and the dashed line is the analytical result. 



X/7t=0.1,X=10, Aj==1.0 
smoothened 
average 




-z=t-v„X 

FIG. 8. The neutral mode contribution to the integrated T>^ 2 (^) at T = for v n = 1, 
A/-7T = 0.1, X = 10, and Ag = 1.0. The thin solid line is for a given realization of disorder, the 
dashed line assumes a measurement with a finite time resolution, and the thick solid line is the 
analytical average. 
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-1 -0.5 _ %=M , x 0.5 1 



FIG. 9. The neutral mode contribution to the integrated Vf 2 (|9|) at T/A ? = 2000 for v n = 1, 
A/-7T = 0.1, X = 10, and = 0.1. The thin solid line is for a given realization of disorder, the 
dashed line assumes a measurement with a finite time resolution, and the thick solid line is the 
result of averaging over 4000 impurity configurations 



There are several requirements which must be met to make the measurement useful. 
First, one must be able to separately (or at least differentially) contact the edges of the 
bilayer system. If each electrode used in the measurement couples identically to both edges, 
one cannot hope to probe the dynamics of the neutral mode. Indeed, the sum of the elements 
of the matrix density-density correlation function for the 331 state ([TTJ) is identical in form 
to that of the Pfaffian state Q2(f). The relative strength of the signal from the neutral mode, 
compared to the charged mode, would be maximized by applying a voltage antisymmetric 
between the layers. Second, the experiment must involve a time-resolved measurement in 
order to distinguish signals which differ by their propagation velocities. Third, the electrodes 
must be close enough together so that the decay of the neutral mode signal with distance 
does not cause it to be undetectable, but they must be far enough apart so that the charged 
and neutral mode signals are well separated in time. 

We believe that the length and time scales needed in a realistic measurement would 
require a careful choice of fabrication techniques. To give more specific estimates, let us take 
the drift velocity of the charge mode to be v c ~ 10 7 cm/s p9| , |20| . In a clean system at a finite 
temperature the strength of the neutral mode signal decays exponentially with distance, see 
Eq. (p3|). The separation between the electrodes here cannot be taken to be much larger than 
the temperature coherence length Lj- ~ hv/T, which is 10~ 2 cm at 10 mK. Assuming the 
neutral mode velocity is smaller than the charged mode velocity by an order of magnitude, 
the arrival time difference between the two modes is on the order of 10 ns. In the presence 
of disorder, the temperature coherence length looses its importance, as one can see from the 
second term in Eq. (|24|). In this case, let us assume a mean free path of £ ~ 1/A^ ~ 10~ 5 cm, 
and require that the electrodes need to be approximately 100£ ~ 10 fim apart for the neutral 
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mode signal to be detectable. The arrival time difference between the signals from the 
neutral and charged modes is then around 1 ns. If an experimental measurement like the 
one described here detected the neutral mode signal, it would conclusively show that the 
v = 1/2 plateau in bilayer systems is the 331 state rather than the Pfaffian. 



III. DETAILS OF THE CALCULATIONS 

In this section we present the detailed calculations of the results summarized in the previ- 
ous section. We begin with a review of the clean edge theory (Section [III A|) , and then discuss 
the addition of a parallel magnetic field (Section [III B|) and disorder (Section [III C| ). Finally 
we discuss the Pfaffian edge (Section |Ti fD|) and the numerical computations performed for 
the 331 edge (Section PTTEf). 

A. Edge Theory of Clean Bilayer Systems 

In this section we review the edge theory of clean bilayer quantum Hall systems with 
interlayer electron tunneling. For a more detailed discussion see I. The edge theory corre- 
sponding to the Halperin state ([!]) contains two chiral Bose fields, Ui(t,x) (i = 1,2), with 
compactification radii Ri = 1 (i.e., u ~ u + 2ir), and equal-time commutation relations: 

[ui(t, x),Uj(t, x')\ = i7iKijSgn(x — x), (30) 

where K is a symmetric, integer-valued matrix which characterizes the topological properties 



of the edge and is completely determined by the exponents in the bulk wavefunction |21 



k= m n ). (si) 

V n m J 

In terms of these fields we can write the charge density and electron creation operators as 
Pi {x) = ±K^d xUj (x), *\(x) = j^e-^\ (32) 



and the Hamiltonian as 

r-L/2 

Ho = I dx 



L/2 



where 



Vij :d x Uid x Uj\ +X (:^i(x)^l(x): +h.c. 



is a symmetric, positive definite (g 2 < v 2 ) matrix which includes the effect of the confining 
potentials and interactions at the edge, and Ao is the interlayer electron tunneling amplitude, 
which we take to be real without loss of generality. The normal ordering is with respect to 
the oscillator modes of the bosonic fields. 



(33) 
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The Hamiltonian and commutation relations can be simplified by the transformation: 



( m 

in terms of which we have 



■y/m — n 



(35) 



!>i(x), 0j(V)] = "i7i"5ijSgn(x — x') 



and 



L/2 
-L/2 



dx 



(36) 



(37) 



where in 0j, the index i runs over the two values i=c, n, which denote the charged and neutral 
modes, respectively, and we have introduced the parameters (3 = \j2{m — n), X = X L~ n , 
the velocities v c>n = {m ± n)(v ± g), and the short distance cutoff a. The Hamiltonian 
separates into a free chiral boson Hamiltonian for the charged mode and a chiral sine-Gordon 
Hamiltonian for the neutral mode. 

For future reference we record the expression for the electron and density operators in 
terms of the newly introduced bosons 



*i,a(aO 



iyj{m+n)/2(t> c (x) =FW>n(z)/2 



L m/2 

, d x (f> c (x) =F — ^d x (p n (x) 

2irJ2(m + n) 2nf3 



(38) 
(39) 



which follow from Eqns. (^) and (|35|). In the remainder of the paper we will suppress the 
subscript on the neutral boson, i.e., <fi = n . 

The time-ordered correlation functions given in Section [11 A| follow from Eqns. (|37|), (p8|), 
and (|39|) ; for details see I. To transform from the time-ordered correlation functions to the 
retarded response functions we note that if the time-ordered function is expressed as 



C(t) = 9(t)C > (t) + 9(-t)C < (t), 



(40) 



where 9(t) is the Heaviside step function, then the corresponding retarded correlation func- 
tion is 



C R (t) = 6(t)[C>(t)-C<(t)}. 



(41) 



B. Parallel Field 



We consider a parallel magnetic field along the y-axis: B = By. This corresponds to 
a vector potential A(^) = Bz~k, where we take the origin of the z-axis midway between 
the layers, whose separation is d. We incorporate this parallel field into our edge theory by 
modifying the charge density operator via the replacement 
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pi(x) i-> pi(x) 



—A x (zi). 



where z\p = ±d/2. Using this along with the definition of the charge density 
transformation (|35|) gives a Hamiltonian 



L/2 
-L/2 



dx 



2tt 



d x 6 + 



47Tt> n 



(42) 
and the 

(43) 



where T = v n Bd(3/2. The second term is a constant, and thus produces only an overall 
shift in the energy spectrum, and we will henceforth ignore it. The inter layer magnetic field 
couples only to the neutral mode, and therefore we will not write the terms involving the 
charged mode explicitly except when considering correlation functions. 

As remarked in Section [[IB , the analysis here also applies to the case where we introduce 
an electric potential difference between the layers instead of an interplane magnetic field. 
A potential difference V e between the layers adds a term V e (pi(x) — P2{x)) ~ V c d x <p to 
the Hamiltonian. This is the same form as the interplane field perturbation (f43|). The 
only difference is that in the case of a potential difference between the layers the density 
operators are not modified as they are in the interplane magnetic field case (R21). However, 



since V^it^x) involves the density fluctuation operator, pi — (p^, see Eq. ([3|), our results 
below for the density- density correlation function apply for either a parallel magnetic field 
or an electric potential difference. 



1. 110 Sequence 



The states in the 110 sequence correspond to f3 2 = 2. In I it is shown that at this value 
of (5 the radius of the neutral boson is R n = 1/ \/2 and therefore we can define a triplet of 
su(2)i Kac-Moody (KM) currents 



J*(x) 

j±rx) 



1 



2ttV2 

j x ± ij y 



dMx) 



2na 



in terms of which the Hamiltonian reads 



H 



B 



L/2 
-L/2 



dx 



27TVr, 



:\3(x)} 2 : +V2TJ z (x) + 2\r 



where we have used the identity 

rL/2 



L/2 



dx :[J z {x)Y: | 7 > 



L l 2 1 2 

dx - : [J(x)] : I7), 

-L/2 3 



(44) 



(45) 



(46) 



valid for any state I7) in the Hilbert space. Next, we can define a new set of currents 
J a (x) = R ab J b (x), which also obey an su{2)i algebra provided R G SO(3). In particular if 
we choose 
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R 



I cos a — sin a 

1 
\ sin a cos a 



(47) 



where sin a = A/yA 2 + r 2 /2, and express the rotated su{2)\ currents in terms of a new 
radius R = l/y/2 chiral boson, 9{x), we then have 



n 



B 



L/2 

-L/2 
L/2 

-L/2 



dx 



dx 



2nv n 



J z (x, 



—v D :(d x 9) 2 : 

47T 



: +2^/A 2 + r 2 /2J z (x) 
1 



7T 



v/2 



a 2 + r 2 /2 a s e 



(48) 
(49) 



Note that this final form of the Hamiltonian is identical to the neutral mode Hamiltonian 
in the absence of a parallel magnetic field with the replacement 



A i-^ A' = a/A 2 + T 2 /2 



(50) 



The above diagonalization of the Hamiltonian also allows us to find correlation functions. 
In particular, consider the two-point function of the density- fluctuation operator, Vij(t,x) 
(0). Using the minimal coupling prescription (|4*2]), the transformation (f4"?D , the definitions 
of the charge density (|39|) , and su(2)i currents (|44[) , we can write 



Pl,2{X) 



1 A A 1 

=o x (p c =F — 
87T 2 (2m - 1) 2tt 



r 



; nX /2 



ig^-^oos(^)} 



Using this expression for the charge density operators in terms of the fields <ft c , 
along with the Hamiltonian (ffHf), we can readily find 



(51) 
and 6, 



W(t,x) 



1 



2(2m - 1) [2tt(x - v c t + ie t )} 2 
1 (1 - a x ) I 



+ 



2(A') 2 [2ir(x - v n t + ie t )}' 



r 2 x2 /2A' x 

h A cos — x 

2 \v n . 



(52) 



The evaluation of the single-electron Green's function is more involved because the elec- 
tron operators (\I/j(x)) cannot be expressed in terms of the fields <ft c and 9. Following the 
method used in I, we can use the independence of the charged and neutral modes to write 
the Green's function, Qij{t,x) (|J), for all states in the 110 sequence as 



Gij{t,x) 



,(110) 



G£\t,x) 



(t, X / 



(53) 



where Gij 1 ^ is the Green's function for the special case m — 1, n — and 



lm-l/2 



,i^/m- ±<t> c (t,x) -i^/m- i</> c (0,0) 



[2tt(x - v c t + iet)] 1 



(54) 
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The decomposition ( |53|) is useful because for the uncorrelated integer 110 state there exists 
a chiral fermion description of the edge theory including tunneling and a parallel field: 



-1/(110) 



L/2 
-L/2 



dx : 



J. 



If we perform the following canonical transformation 

- J<P<7v ( ^+ 



^1 

^2 

then the Hamiltonian ( j55f) becomes 

-L/2 



A 



where sin(2<^) = — — 
A 



(55) 



(56) 



L/2 



<ix : 



(57) 



By bosonizing according to ip±(x) = e 1 ^^ / y/2na, and defining 9 + = (0+ + 0_)/a/2, and 
= (0+-0_)/ v / 2-v / 2A / xK, we can exactly evaluate the single-electron Green's function 



g (110) (t,x) 



1 



2ttJ (x — v c t + ie t )(x — v n t + ie t ) 



exp 



r 



x 



i -^o z + Xa x — 



V2 



(58) 



Combining Eqns. (|53|), (|54|) , and (|58|) , we finally arrive at 

1 



£(M) 



[27r(a; - v c t + ie t )] m '^sj2n(x - v n t + ie t ) 



x 



i ( r 



11 cos(A'a;/f n ) + — + Ac*) sin(A / x/'u n ) 



(59) 



2. 331 Sequence 

The states in the 331 sequence correspond to f3 2 = 4. In I it is shown that at this value 
of f3 we can fermionize the neutral boson using 



J4> 



J2(j> 



ipd x ijj. 



r ' 2vr ^ ' r r '' 2vra 2 
With these identities, the neutral mode part of the Hamiltonian ( fi3| ) becomes 



(60) 



n 



B 



L/2 

-L/2 
L/2 

-L/2 



dx : 



cfe : 



-iv n ^d x ij; - i^-{ip ] d x ^) ] + ipd x ip) + T^V 

Z7T 



■^lxAxi - ^2X2^x2 + ^xiX2 



(61) 
(62) 
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where we have written the chiral Dirac fermion in terms of its Majorana components: ip(x) = 
\Xi{ x ) + ^X2(a ; )]/v / 2? where %■ — Xi> an d recalled t>i )2 — v n ±X/ir. The tunneling term splits 
the velocities of the two Majorana fermions and the parallel field term couples them. 

The Hamiltonian is quadratic and hence readily diagonalizable. If we take antiperiodic 
boundary conditions for the Fermi field, ip(x + L) — —ip(x), and expand in Fourier modes 
according to 



where k G (27r/L)(Z + 1/2), the Hamiltonian (|6T| ) becomes 



Ti 



B 



E : 



v n kc\c k + Tc\c k + ^-{c\c ] _ k - c k c- k ) 

Z7T 



E : ( 4 c -k 

k>0 



v Q k + r Xk/7T 
Xk/rr v n k — T 



Cfe 

c -k 



Employing the Bogoliubov transformation 



(63) 



(64) 



c -k 



cos a k — sm a k 
sinafc cosafc 



where a_ 



— afc, the Hamiltonian is diagonalized via the choice 

Xk 



tan(2afc) 



(65) 



(66) 



along with the restriction 2a k G [— ir/2, vr/2], required to produce the correct spectrum in 
the limit A — > 0. This yields 



U B = E [v*k + sgn(r) v /(AA;/7r) 2 + T 2 
= J2e(k) :a\a k : . 



: a k a k : 



(67) 



To calculate correlation functions, we first use the transformation flBBD , along with the 
expressions for a k fl66| ) and V'O 3 ') (p3| ) to express the Fermi field in terms of the mode 
operators which diagonalize the Hamiltonian 



..ikx 



ip(x) 



1 & 

W\ [(AA;/7r) 2 + r2]i/4 



X 



(AA;/vr) 2 + T 2 + \T\ a k -sgn(TXk)\ J (Xk/n) 2 + T 2 - \T\ a 



To compute the single-electron Green's function (0) we use the transformation 
fermionization fl60|) to write the electron operators (j32|) as 



1 



(2vra)( m - 1 )/ 2 



(68) 
']) and the 



(69) 
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Since the charged and neutral modes are not coupled, we find 



^ ^(x-J+^r-' ^' (70) 

where we have defined the matrix 

c a r \ = ( s )^(o> o)) (mt, x)m o)) \ 

The correlation function of the fermionized neutral mode can be reduced to quadrature. 
Using Eqns. ( |STD and ( |B"8| ) we obtain (in the limit L — > oo) 

_ -< / i + (i/t)(0/ar) {\/irTt){&/dx8r) \ 
4?r {(x/nrt^dP/dxar) i-{i/t){d/dv) J > 1 J 

where there remains the integral 

I = J dk [sgn(t) cos(fcr — e(A;)t) + zsgn(fc — k F ) sin(/cx — e(A;)t)] , (73) 

and the "Fermi momentum", defined by e(kp) = 0, is given by kp = —V/ ^Jviv 2 - To simplify 
the result, we first change the variable of integration to uo = e(k) for k > kp and uj = —e(k) 
for k < kp. With these substitutions and some algebra, the matrix Green's function can be 
written 



— 2 

QJt,x) = sgn(t) 

VlV 2 



v n TL - -(A P x {t,X) + sgn(t) (-1 - v n <A -P T {r,X) 



-n;it; 2 rVP(T,X) 

where we have defined X = x/viv 2 , r = sgn(t)(t> n X — t), and the function 



(74) 



^,^r^^, (75) 



where f«(u;) = yviv 2 T 2 + (X/tt) 2 uj 2 , and the subscripts on P in Eq. (|T3| ) denote partial 
differentiation. Although we have been unable to evaluate P(t, X) explicitly, the real part 
of this function can be calculated in closed form. This is discussed below in Section [III C| 
when we consider the retarded version of Q^. 

Turning now to the density-density correlation function, we can use the fermionization 
(BO) to write the density operators (p9]) as 



p 1>2 = - T d x <p c T \ : 1> V : • (76) 

4:7T\/m — 1 2 



Since the Hamiltonian in terms of ip ( |6T|) is quadratic, Wick's theorem holds for this field 
and the density two-point function can be expressed in terms of the single-particle Green's 
function. Using ( [76]) and @ one finds 
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W(t,x) 



1 



4(m - 1) [2tt(x - v c t + ie t )] : 



+ - — - — '-detg^x). 



(77) 



Along with Eq. ([74]), this gives an expression for the density-density correlation function in 
terms of the single function P(t,X) and its derivatives: 



iV(t,x) 



'II - a 3 



4(m - 1) [27r(x - v c t + ie t )} 2 Av x v 2 



P 



v- ^p r 2 + ^ 2 r 2 p 2 



(78) 



C. Disorder 

We now consider the effects of adding several types of disorder to the Hamiltonian of 



the bilayer system (33). We begin by considering the relevancy of various random terms 
within a renormalization group (RG) analysis. We then present exact results for the 110 
and 331 sequences, concentrating on the retarded density response function because of its 
relevancy for experiments. The 110 case is solved by using an SU(2) gauge transforma- 
tion to separate the Green's functions into products of clean Green's functions and terms 
involving only the random fields. After this step, the disordered problem is shown to be 
equivalent to a spin-1/2 particle in a random magnetic field and the disorder averaging is 
performed non-perturbatively. The solution for the 331 sequence involves an exact summa- 
tion of the disorder-averaged perturbation theory which is possible because of the chirality of 
the system. In Appendix |A| we present an alternate method for obtaining disorder-averaged 
correlation functions for the 331 sequence based on the spin analogy. 
If we consider a general perturbation to the Hamiltonian of the form 

f L/2 

U 5 = / dx((x)0(x), (79) 

J-L/2 

where 0{x) is an operator of scaling dimension 6 and ((x) is a Gaussian random variable 
with variance A, i.e., ({x)({x') = A5(x — x'), where the bar denotes disorder averaging, 
then a lowest-order perturbative RG analysis gives p2 



^ = (3 - 26) A, (80) 

where the short distance cutoff increases as I increases. 

Consider first the possibility of disorder in the velocity-interaction matrix V in the Hamil- 
tonian (|33|). Since the V matrix multiplies an operator of scaling dimension 5 = 2, we see 
from the flow equation flSOl) that delta-correlated disorder is RG irrelevant. Therefore we 
can ignore randomness in the V matrix and interpret the values appearing in Eq. ( |34D as 
disorder-averaged mean values. Note that the V matrix must be symmetric, since it multi- 
plies a symmetric operator in the Hamiltonian, and the positive-definiteness of V is required 
for the Hamiltonian to be bounded from below. However, the assumption that Vn = V22 is 
made for technical reasons and it can be relaxed to V\\ = V22, a weaker criterion. 

Next we turn to the case of disorder in the tunneling amplitude A in the Hamiltonian 
(|TF|). The scaling dimension of the tunneling operator it multiplies is 5 = (3 2 /2. Using 
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Eq. (|80D , we find that disorder in the tunneling amplitude is relevant for the 110 sequence 
(i.e., p 2 = 2), and irrelevant for the 331 sequence (i.e., p 1 = 4). The former case will be 
discussed later in this section and the latter case in Appendix [B]. 

Finally we consider adding random scalar potentials which couple to the edge charge 
densities in each layer: 



-L/2 rL/2 

dx [£x(x)px(x) +£ 2 (x)p2(x)] = 

-L/2 J-L/2 



^c(x)^d x (f) c (x) + i n (x)^-d x (t)(x) 



(81) 

where we have used the definition of the charge density operators (|39"D, and 



i = i , 1 , (6 + 60, ^ = i(6-eo- (82) 

'2(m + n) P 



If we assume that Ci,2( x ) are independent Gaussian random variables then so are £, CjTL (x). 
Since these terms involve operators of scaling dimension 5 = 1 they are relevant perturba- 
tions for both the 110 and 331 sequences. 



1. 110 Sequence 



In the discussion above we found that for the 110 sequence (/3 2 = 2), disorder both 
in the tunneling and in the scalar potential terms is relevant. We therefore consider the 
Hamiltonian 



Ti 



D 



L/2 
-L/2 



dx 



£ :(<9,0 C ) 2 : +^ -.(dAf: +-L (\(x)e^ + A*(x)e~^)) 
47T 47T Ana v ' 



1 1 

+ (cWrUW + ^n(x)— 9^0(0;) 

Z7T Z7T 



(83) 



where \(x) is a complex random tunneling amplitude. The presence of disorder breaks 
translation invariance and hence the current algebra method used to solve the clean problem 
in Section 003 is of no use because the transformation J a (x) = R ab J b (x) must be a global 



rotation to map between sets of KM generators. However, an alternate approach to the 
problem developed in I is useful in the disordered case. We add to the Hamiltonian ( |8"3"D an 
auxiliary free chiral boson ((f)) with a velocity equal to the velocity of 0: 



rL/2 1 

Hd | — > 7~(-d + / dx — v n :(d. 
J-l/2 An 



2 . 



perform the canonical transformation 



J_/l IWfc 

v^u -i)\e 2 



(84) 



(85) 



and then fermionize according to 
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(86) 



The details required to make this mapping rigorous (i.e., compactification radii, topolog- 
ical charges, Klein factors) are discussed in I. The result of this procedure is a quadratic 
Hamiltonian 



H D = f L/2 dx \^-v c :(6W> C ) 2 : +Ux)^d x <p c (x)+ : (-iv u tfd x * + v n B a (x)¥a a V 

J-L/2 L47T Z7T v 



(87) 



where we have defined 



(^j) , B(x) = ±(-te[\(x)),-^[\(x))M*)/yfy, (88) 

and the index a runs over x, y, z. The fermionic part of this Hamiltonian describes a pseudo- 
spin-1/2 fermion coupled to a random SU(2) gauge field. 

We now perform a change of variables that absorbs the disordered terms into the defini- 
tions of the field operators. For the charged mode we define 

r)(x) = <p c (x) + - / dyUv), (89) 



and for the neutral mode we use an SU(2) gauge transformation 

V(x) = S(x)V(x), (90) 
where S(x) G SU(2) is a solution of the matrix differential equation 

= -iB a {x)a a S{x). (91) 
With these definitions, the Hamiltonian (B7T) becomes 



d = f L/2 dx [-U c :(d x 7]) 2 : -J—g-iv n :¥d x V: 

J-L/2 L47T AlTVc 



1 ,n 1 



(92) 



Since the second term only involves the disordered scalar potential, it does not affect corre- 
lation functions and will henceforth be neglected. Note that in going from Hamiltonian ( |87| ) 
to Hamiltonian (|92D we have used a gauge transformation on a chiral Fermi field (0) with- 
out accounting for the chiral anomaly. This is valid because the gauge field is a quenched 
random variable; in this case the anomaly associated with the chiral gauge transformation 
(p0|) cancels in the average. 

Our primary goal is to understand the behavior of the density-density correlation function 
in a given sample, i. e., for a given realization of disorder. We begin by expressing the density 
operators (50) in terms of the fields r\ and with the help of Eqns. (B3F), (BBT) , and (B9I): 



p h2 (x) = , 1 (d x r]{x) - -&xj) =F - : ¥{x)a z ^{x) : . (93) 



2nJ2(2m- 1) 
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Using this expression in the definition of the density two-point function (|D we have 

1 (H_+a x ) (1 — a x ) 

- iV(t, x, x') = - rj—, — — - tttt - - -tr (a z G(t, x, x')a z G(-t, x', x)) 

v ' ' ; 2(2m-l) [2n{x-x' -v c t + ie t )} 2 4 v v ' ' ' v ' ' ;; 

(94) 

where we have used the single-particle matrix Green's function 

G i:i {t, x, x') = -z(T^(t, x)^J(0, x')). (95) 

We have explicitly included two spatial arguments in these correlation functions because 
of the lack of translation invariance in a given realization of disorder. Note that although 
the charged mode, C , is coupled in ([81]) to a disorder potential, £ c , the charged-mode part 
of the above correlation function fl94|) is identical to the result in the absence of disorder. 
This result is true for every realization of disorder and is essentially equivalent to the loop- 
cancelation theorem which states that for linearly dispersing fermions (e(k) cx k) in 1+1 
dimensions the connected n-point function of the density operator vanishes identically for 

n> 2 jag. 

To determine the effect of disorder on the neutral mode we first note that the differential 
equation (|91|) has a solution in terms of a coordinate-ordered exponential 

S(x) = T y exp (-i f dyB a {y)a a ^ , (96) 

where T y is the y-ordering operator. Since the matrix S(x) can be taken outside quantum 
expectation values, we can express the Green's function of the ^ field ( j95|) in terms of the 
Green's function of the (free) ^ field and thus write 



G^x.x') = — ^—Uifax!), (97) 

27r(x — x' — v n t + ie t ) 



where 



U(x, x') = S(x)S j (x') = T y exp (-i £ dy B a (y)a a ^j (98) 

is a unitary matrix. Using Eq. ( p7|) we find that the neutral mode part of the density- density 
correlation function is proportional to 

-tr (a z G(t, x, x')a z G(-t, x' , x))- — -rtr (a z U(x, x')a z UUx, x')) , (99) 

[2tc(x — x — v n t + ie t )\ 2 v ' 

where we have used the property U'(x, x') = U(x', x), which follows from the definition (|9*8"D . 
In Eq. (p9|) we have written a correlation function in the disordered system as the product 
of the corresponding function in the clean system and a factor which depends only on the 
random potential. 

If we interpret the coordinate y appearing in the definition of U as a fictitious time, then 
the matrix U (x, x') is exactly the time evolution operator between times x and x' for a zero- 
dimensional system with time-dependent Hamiltonian B a (y)a a . This is the Hamiltonian for 
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a spin-1/2 object in a random magnetic field B a (y). The quantity appearing in the trace in 
Eq. ( p9|) can then be interpreted as the (S z (x)S z (x')) correlation function for this spin. 

To understand the behavior of the density-density correlation function in a given sample 
we will first calculate its disorder-average, which involves averaging the quantity appearing 
on the r.h.s. of Eq. fl99|). Toward this end, consider the following vector quantity 

(100) 



F a (x; x') = tr (c/t(z, x')a a U(x, x')a z ) 
By differentiating with respect to x we find that this is a solution of the differential equation 

dF a (x-x') 



dx 



M ab (x)F b (x-x') 



(101) 



where M ab (x) = —2e abc B c (x), subject to the boundary condition F a (x';x') = ti{a a a z 
2S az . The solution of this differential equation can also be written as 



F a (x:x') 



T y exp (J ' ^ dyM(y) 



ab 



F b (x'-x'). 



(102) 



We have expressed the quantity we desire, F z (x; x'), in terms of a single coordinate-ordered 
matrix exponential, whose disorder average we now show can be readily evaluated. 

We assume the tunneling amplitude and scalar potential are delta-correlated Gaussian 
random variables and denote the mean and variance of B a (x) by fi a and A a , respectively. 
Since the exponential appearing in Eq. (|102p is ordered in y, and the elements of M(y) are 
independently distributed for each y, we can consider breaking up the interval [x', x] into N 
intervals of length e = \x — x'\/N and then taking the limit N —>■ oo. Therefore we can write 



T y exp( [ X dyM(y))= lim \ f dB P(B)e ssn{x - x ' )eM 

\Jx' J N-*oo J 



N 



(103) 



where the probability distribution is 
P(B)rfB 



exp 



1 3 e 

a=l 



dB. 



\\ (2vr) 3 A a: A^A 2 

Expanding the e B sa-( x - x ') eM factor in Eq. ( |103| ) and performing the integration gives 

2\x — x'\ 



(104) 



T y exp 
where 

W = 



dyM{y) 



lim 

N^oo 



1 - 



N 



-W + O 



1 

iV2 



A' 



exp(-2\x-x'\W) } (105) 



A y + A z 
A z + A X 
A x + A y 





( o 


// 2 




+ sgn(x — x') 


-fi z 











—/jL x 






(106) 



From Eqns. ( |100| ), ( |102| ), and ( |105| ), we finally arrive at 

ti(W(x,x')a z U(x,x')a z ) - <> \~- 2 \*-> r ' n 



(107) 
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While it is in principle possible to evaluate the exponential of W for arbitrary /i a and 
A a , for simplicity we shall restrict ourselves to the case of a real tunneling amplitude. If the 
tunneling amplitude has mean A and variance A^, and the disordered scalar potential has 
mean zero and variance A^, then from the definition of B(x) ( jS8f) we have 



A 

V = , 



A a 



Aa 



fi y = 0, // 
A y = 0, A 2 





A? 

2vf 



In this case 



W 



a A A + A e /4 + g . A 



4f2 



ia" sgnix — xl — 



(108) 



(109) 



and from Eqns. ©, @, and fllQ7|) we find 

1 (H. + o*) 



iV(t, x, x') 



2(2m - 1) [2vr(x - x' - v c t + ie t )} 2 

(1L-a x ) ( 2\x-x'\ 



1 

+ ^ 



2 [2ir(x - x' - v n t + ie t )]' 



exp 



x 



cos 



2\x — x'\X 



A, 



'2\x-x'\\ s 



(Aa + A € /4) 



4f n A 



^ sm 



(110) 



where 



A = 



\ 



A 2 



V 2 

4v„ 



(111) 

and fl41|) we arrive at 



Transforming to the corresponding retarded function using Eqns. 
Eq. @. 

The question remains as to what the behavior of T>(t, x, x') is in a given sample. Instead 
of attempting to evaluate higher moments of this correlation function, we shall exploit the 
fact that it can be expressed in terms of the single-particle Green's function, whose second 
moment we will compute. Writing out the trace in Eq. (|94]) explicitly we find 



tr (a z G(t, x, x')a z G(—t, x', x)) = Gu(t, x, x')Gn(—t, x , x)—Gi 2 (t, x, x / )G 2 i(— t, x , x) 

-G2i(t, x, x')G 12 (-t, x', x)+G 2 2(t, x, x')Gzi(—t, x', x). (112) 



From Eq. 



and the result 



U (x, x') 



exp 



\x — X 
2v 2 

n 



(A A + A e /2) 



11 cos — (x — X 



ia x sin 



A 



x — x) 



(113) 



which can be obtained via the same procedure used to evaluate the average in Eq. ( |103| ), we 
find the disorder-averaged single-particle Green's functions are also exponentially decaying 
in space 



29 



G\i(t, x, x') = G22(t,X,x') 



G 12 (t, x, x 1 ) = G 2 i(t, x, x') = i 



r COS (x — x') 

2n(x - x' - v n t + ie t ) \v n y 

e -\x-x'\(A x +A ( /2)/2v2 / A 

— r sin — (x — x') 

2n(x — x' — v n t + itt) \v n 



(114) 
(115) 



To investigate whether this exponential decay is an artifact of the disorder averaging (i.e., it 
arises from averaging over random phases) , or whether we expect it to hold in a given realiza- 
tion of disorder, we compute \Gij(t, x, x')\ 2 , which is clearly insensitive to phase fluctuations. 
From Eq. (^) we see that 



\Gij(t,x,x')\< 



\2tt(x — x' — v n t + ie t )\ 2 
By the unitarity of U(x, x') we have 

u n uli + u 12 u 2 \ = i, 



(x, x'). 



(116) 



(117) 



where we have suppressed the spatial arguments. Using the explicit form of U that follows 
from Eqns. (0), (0), and 3m[A(x)] = 0, one may show a v U T a y = U\ which implies 



u^uli = u 22 ul 2 , U 12 W 21 = U 21 U-_ 



12- 



Finally, from Eq. ( |107| ) we have 



U u Uli - U 21 Ul 2 - U 12 U j 21 + U 22 U\ 2 = 2 



-2\x-x'\W 



118) 



(119) 



Eqns. (|117|) -( p!T9|) are four equations in four unknown quantities which can be solved to yield 



UnUl = U 22 U\ 2 = - (l 



-2\x-x'\W 



U\ 2 U 2 i — U 2 \U 12 



-2\x-x'\W 



(120) 



This result, combined with Eq. ( 116 ) gives the disorder- averaged absolute- magnitudes of the 
elements of the single-particle Green's function 



I Gij (t, X , x ) I 



1 



2|27r(x — x' — v n t + iet) 



2\x — x'\ y 



cos 



+ 



A 



1 + (~l) i+j exp y 
- . ( 2\x — x'\\ 



2\x — x'\ 



(A A + Af/4) 



Av n X 



— sm 



;i2i) 



The structure of this result is interesting. Each \Gij\ 2 is the sum of two terms, one of 
which is identical to the square of the free (i.e., \(x) = (, n (x) = 0) Green's function, and 
the other of which has spatial oscillations at the shifted frequency and an exponential decay 
in space from the disorder. The fact that |Gy| 2 has a long-ranged part (i.e., a term that 
decays algebraically rather than exponentially) indicates that in a given sample G^ is not 
exponentially damped, and thus from Eq. ( p.!2| ) we can conclude that the neutral mode 



30 



portion of the density-density correlation function is also long-ranged for a given realization 
of disorder. We therefore expect that in a given sample the density two-point function has 
the structure of the disorder-averaged quantity ( |110| ), without the exponential decay in space 
of the neutral mode piece. As a function of Ax = x — x' we expect two peaks at the points 
Ax = f C)n t, with the second peak modulated in space at a frequency that varies with the 
local scalar potential [see Eq. (11111) 1 ■ 

2. 331 Sequence 

In the analysis at the beginning of this section we determined that for the 331 sequence 
(f3 2 = 4), only disorder in the scalar potential terms was a non-irrelevant perturbation. We 
therefore consider the Hamiltonian 



fL/2 

Hd = dx 

J-L/2 



+ £c{x)—d x (f> c (x) + i n {x)—d x (j){x) 

Lix 111 



(122) 

which, with the help of the fermionization ( j6Qj ) for the neutral boson and the transformation 
( jS9|) for the charged boson, can be written 



fL/2 

Hd = dx 

J-L/2 



-^-v c (d x rj) 2 - -vxxAxi - ~v 2 X2d x X2 + i£n(x)xiXz 



(123) 



where we have dropped the constant ^ term. The charged mode portion of this Hamiltonian 
is identical to that of Q52] ) for the 110 sequence, and hence all results pertaining to the charged 
mode can be imported from the previous discussion. The Hamiltonian is now quadratic, 
however the lack of translation invariance prevents us from employing the method used in 



Section |III B| , where we essentially solved the special case in which £ c (x) is independent of 
x. In addition we cannot absorb the disorder into the definition of the field operators via a 
gauge transformation as we did for the 110 sequence, because the Majorana fields are real 
and therefore neutral. However, we can still separate the disorder and quantum expectation 
values by explicitly constructing solutions to the Heisenberg equations of motion. This 
procedure involves some technical subtleties not present in the 110 solution, and is discussed 
in detail in Appendix [A[ In this section we use a different approach; we find the disorder- 
averaged correlation functions of the above Hamiltonian by an exact summation of the 
disorder-averaged perturbation theory. 

The chirality of the fermions in the neutral mode part of the Hamiltonian allows a great 
simplification in the structure of the diagrammatic perturbation theory in powers of the 
disorder potential £ n {x)- This was first noted by Chalker and Sondhi in the context of a 
single-particle description of the edge ||15|| . Consider the matrix Green's function of the 
Majorana fields for the free case, i.e., £ n ( x ) — 0: 

g§\t,x,x') = -i(Txi(t,x)xj(0,x')) = — ^ , (124) 

J 2ir{x — x — Vjt + te t ) 

Fourier transforming with respect to time gives 
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Mr, , ~ „'\ — x I A+J^t 1 



glj'(uj,x,x') = 5ij j die 



2ti{x — x' — Vjt + ie t ) 
5 ..l e M*-x')/v 3 \e{-uj)e(x' -x)- 6{uj)6{x - x')\ . (125) 



To obtain this result note that the integrand has a pole in the complex t plane at t — 
(x — x')/vj + 0(e). Therefore for x — x' > 0, 3?et > at the pole and hence the pole lies in 
the upper half plane while for x — x 1 < 0, 3?e t < at the pole and it is therefore in the lower 
half plane. We find that for positive frequencies (u> > 0) the function vanishes for x — x' < 0, 
while for negative frequencies (cu < 0) it vanishes for x — x' > 0. 

Next consider the single-particle Majorana Green's function with the disorder potential 
present in Eq. ( |123| ). Working perturbatively in powers of £ n (a0 we have 

oo . 

g n (uj,x,x') = Y. d yi--- d y2ngn\^,x,y 1 )^ n (y 1 )g^(u,y 1 ,y 2 )^n(y2)---gn\uJ,y2n,x'), 

n=0 J 

oo , 

g u (uj,x,x') =iY. d Vi- ■ ■dy2n+ign\uJ,x,y 1 )£ ri (y 1 ) . . .g$(u,y 2n +i,x'), (126) 

n=CT 

with similar expressions for the remaining components. Each time the particle scatters off 
the impurity potential its velocity changes from v± to v 2 or vice versa. When we disorder 
average the above equations, we must tie together insertions of £ n (x) in all possible ways for 
each term in the sum. We first observe that because the scalar potential is time- independent, 
io is conserved and hence it has the same sign in every propagator. This, along with the 
chirality of g(°\tu, x, x') evident in ( |125[) implies that any disorder-averaged diagram in which 



impurity lines cross vanishes identically. Therefore, for each term in the sum in (|126|) there 
is a single nonzero disorder-averaged diagram, i.e., the one in which successive insertions of 
£ n are pairwise contracted. The resulting series can be summed to give 

g i3 {uo, x, x') = 8..L e M*-*')/v j(r As\ x - x >\/2v 1 vi jp^ytf _ x ) _ e(uj)6(x - x')} . (127) 

We see that the disorder-averaged Green's function retains the chiral structure of the free 
Green's function. While the exponential decay of the function involves the geometric mean 
of the two Majorana velocities, y / viv 2 , the frequency dependence only contains Vi@) for 
#ii(22)- This is a direct consequence of the fact that all terms with crossed impurity lines 
vanish when we disorder-average the chiral Green's function. Therefore, for the averaged 
single-particle Green's function, if the particle begins propagating with velocity v± it never 
propagates with velocity v 2 , the other velocity enters only through the density of states when 
scattering off the potential. Transforming back to the time domain we find 



e -A ( \x-x'\/2viv 2 

27i(x — x' — vd + ie t ) 



gJt, x, x) = 5 ij7;z7Z — - — — (128) 



From the relation between the Dirac and Majorana fields, ij) = (xi + ^X2)/v / 2, we find 
that the Green's function of the neutral mode fermion, (|7lD , can be obtained from g via 
a unitary transformation 



32 



g f (t, x, x') = Og{t, x, x')0\ where O = -j= Q % _ J . (129) 



Thus from Eqns. (|70|), ( 127 ), and ( 129 ) we find, after Fourier transforming, that the single- 
electron Green's function for the 331 sequence in the presence of disorder is: 



11 

G[t,x,0) = —- — , -— r - 

[2tt(x - v c t + ie t )\ 2 



1 + a x 1 - a x 

+ 



2tt(x — Vit + ie t ) 2ir(x — v 2 t + ie t ) 



-At.\x\/2 Vl V 2 



(130) 



We see that for the single-electron Green's function the velocity split of the neutral mode 
remains in the presence of disorder, but the function acquires an exponential decay with 
distance. 

We next consider the calculation of the density-density correlation function. As a first 
step toward understanding the behavior of this correlation function in a given sample, we 
will calculate its disorder-average. We can use the transformation ( [59] ) and the expression 
for the Fermi density operator in terms of the Majorana fields, :^ r i\):=: 1X1X2 to write the 
density operators (|7^) as 

Pi,2 = - — I- = (d x r) - — £ c ) =f ~ -XiX2-, (131) 

From which we find 

- &(t, *, *0 = 77 1 1U9 7 ( * + g ? ■ • ™ - ^-r^- D & *> ( 132 ) 

4(m — 1) [27c{x — x — v c t + tet)\ 4 
where we have denoted the neutral mode contribution by 

— iD(t, x, x') = (T :xi(t,x)x2(t,x)::xi(^,x')x2(0,x'):). (133) 

To evaluate D, we first Fourier transform with respect to time in order to exploit the chi- 
rality in the mixed frequency-space domain. Computing the expectation value in Eq. ( |133| ) 
using Wick's theorem and taking the disorder-average then gives: 

- W(u, x, x') = -iJdt e wt D(t, x, x') = J ^-F (^^, ^7^, x, x'^j , (134) 
where we have defined 



F(ui, uj 2 , x, x) = guiux, x, x')g 22 {-u 2 , x, x') - g 12 {uJi, x, x')g 21 (-u 2 , x, x'). (135) 

If one uses Eq. ( |126| ) to write the single-particle Green's functions in the above expression in 
terms of the free Green's function, g^\ and the disorder potential, £ n , one finds upon disorder 
averaging that the chirality of Eq. ( |125| ) implies that all non- vanishing diagrams are of the 
form of ladder diagrams with the legs of the ladder constructed out of the disorder-averaged 
propagators g^. 

The legs of the ladder are given by 
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h(ut,u 2 ,k) = J dxe tkx g u (uj 1 ,x,0)g 22 (-uj 2 ,x,0) 

e{ux)e{-u 2 ) e{-u x )e{u 2 ) 



ViV 2 



k - u\jv\ + io 2 /v 2 - iA^/viv 2 k - U\/v\ + uj 2 /v 2 + iA^/viv 2 



(136) 



which was evaluated with the help of equation ( |127| ) for ~g^. The segments of the ladder 
shown in Fig. alternate between h(ui,u 2 , k) and h(—u 2 , —Ui, k). Performing the ladder 
sum gives 



F(ui,u 2 ,x,x') 



dk e ik( x -x') 
2tt 



h(u>i, a> 2 , k) - Ag/t(o>i,a> 2 , k)h{-uj 2 , — o>i, fc) 
1 — A|/i(o;i, u 2 , k)h(—u 2 , —UJ\, k) 



(137) 



F((a l ,G} 2 ,x,x') = 



CO, 



-co 2 



00, 



CO, 



^ — ^ 



i-C0„ 



1 -CO- 1 

I I 



-> — 



v, 



y 5 



FIG. 10. Ladder sum for the neutral mode density-density correlation function. The thin solid 
lines represent g n , the thick solid lines represent g 22 and the dashed lines represent the disorder 
potential and carry a factor of Ag. 



Using the result for h(uJi,uj 2 ,k) in this expression for F(ui,u 2 ,x,x') gives, after some 
algebra, 



F{ux,u 2 ,x,x') 



dk „• 



ik(x— x') 



V\V 2 J 2tt 



0(u;i)0(-u; 2 ) 



k-w(A^) 



[k - z + (As)][k - z-(A c )\ 
k - w(-A^) 



[k - z+(-At)][k - z-(-Ae)) 

where we have defined the parameters 

w(A^) = ui/v 2 - u 2 /vi + 2iA i /v 1 v 2 , 
1 



*±(A 5 ) = 



VlV 2 



V l+ V 2( \ . -A _L ■ A2 (^1-^2)% s 2 

— - — (wi - w 2 ) + *A C ± iW A| (wi + ^ 2 ) 2 



(138) 

(139) 
(140) 



From the expression for z±(A^) we see that in the first term in Eq. (p. 38 ) both poles are 
in the upper half plane while in the second term both poles are in the lower half plane. 
Performing the ^-integration by the residue theorem and using the resulting expression for 
F in Eq. ( PI gives: 



iD(u, x, x') = [6(— lo)6(x — x) — 6{u)6{x — x')\e 



iv n LoX-A^\X\ 



V\V 2 



duo' 
4tt 



A, 



cosh(A(fj')X) + -=^-smh(A(u')\X\) 



(141) 
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where we have used the rescaled coordinate X = (x — x')/viV2, and defined 



A(w) = 



'Af-^i-^ 2 



V^ 2 



(142) 



The experimentally measurable quantity is the retarded density-density correlation func- 
tion. Since Eq. ( |141| ) is not in the time domain, we cannot use Eqns. (flCf ) and fl4l|) to obtain 
D R (t,x,x') directly. However, in uj-k space we have the relations 



D R (u, k) = D(uj, k), 3m D R (u, k) = sgn(cj) 3m D(u, k), 
which in turn yield 

W(u,x) = e(w)D(w,x) + 6(-u)\p(w,-x)]*, 



(143) 



(144) 



where the asterisk denotes complex conjugation. Using Eqns. (|141 ) and ( 144 ) and Fourier 
transforming back to the time domain we find: 



D R (t,x,x') 



9(x — x') 

00 duj 



2vr 



-AtX 



v„X-t 



i'jj(vu X -t) 



cosh(A(o;)X) + 



A ^ sinh(A(cj)X) 



A(w 



(145) 



where V denotes the principal value. The remaining integral is computed in Appendix 0, 
and from Eqns. ( |C7| ) and ( |C8| ) we find 



Z7T 

V 



-A£x/v-lV2 



(vi - w 2 )a; 
\ / x/viv 2 



{vi5(x - vit) - v 2 5(x - v 2 t)} + 



2 Vl v 2 v 1 



v n x/viv 2 - t 



h 



A/tt 



71 T 



V 



TT 



(146) 



where z = (t — x/vi)(x/v 2 — t), and are Bessel functions of imaginary argument. Note 
that as expected this function is real. This result can also be found using the formalism in 
Appendix |A|. 

We have computed D R , but what we are really interested in is the behavior of D R in a 
given realization of disorder. A surprising feature of D R is that it has a (principal value) 
singularity at the mean arrival time to — v n x/v\v 2 . An immediate question is whether this 
singularity is present in each sample, and if not, how does it arise in the average. 

To investigate the behavior of D R in a given realization of disorder we have adopted 
several approaches. First, as in our analysis of the 110 sequence, we use the fact that D R 
can be expressed in terms of single-particle Green's functions, whose second moments we 
evaluate. Second, we show that in a given sample D R (t,x) exhibits an exact symmetry 
about the point t = to- Finally we shall consider the behavior of the correlation functions 
in some simple model potentials. 

The relation between the time-ordered correlation functions D and gij follows from the 
definition ( |133| ): 
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D(t, x, x) = idet g(t, x, x ). 



(147) 



Using Eq. (|40|) , we then conclude from the above equation that D > = idet g> and D< 
idetg < . Thus from Eq. ( f41~|) we have 



D (t, x, x) =i6{t)[detg > {t,x 1 x) - det g < (t, x, x')]. 



(148) 



If we define the advanced correlation functions g A = 9(—t)(g < — g > ), then one can show 
g < = g — g R and g > = g — g A , which when substituted into Eq. ( p.48]) give 



D R = i9(t)[(g u gg-g 12 g£_) 

+ {9ll922 - gi292l) - {.911922 - 9l292l)}- 



(149) 



We see that the retarded density-density correlation function in a given sample can be 
expressed in terms of the time-ordered and retarded single-particle Green's functions. 

We have already evaluated the average time-ordered single-particle Green's function, see 
Eq. ( |128j ). Using Eqns. ( |4"0"| ) and ( |4"TD we find for the corresponding retarded function: 

(150) 



g%{t,x,x') = -i8ij6{t)8{x - x' - Vi t)e- A ^ x - x ' )/2viV2 . 



Next consider the absolute squares of the single-particle Green's functions. The disor- 
der average of gij(t,x,x')g*j(t' ,x,x') can be evaluated by the same diagrammatic proce- 
dure used to obtain the average density-density correlation function. The disorder average 
of g R (t,x,x')g R *(t',x,x') can be readily computed using the formalism presented in Ap- 
pendix [A]. 

Omitting the details of these calculations we find: 

1 e" A «l x l /■«> du , 



g ij (t,x,0)glJt',x,0) 



% cosh(AHX) + 



2tt e + iAt 



D iu}(v n X-t)-(e+iAt)\uj\/2 



27T 



sinh(A(cu)X) 



(151) 



where t = (t + 1' )/2 is the average time, At = t — t' is the time difference, and e > is an 
infinitesimal regulator. For the retarded functions we find 



gft(t,x,0)gR{t',x, 0) 



6{At)e 



-A e \X\ 



S(x/vi —t) + 



g&(t,x,0)g&(!f,x,0) = ?®5(At)e-W X 



VlV 2 



2A 



2A 


6{z 


»■( 


A, 

A/tt 
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x/v 2 — t 
\ t — x/vx 



(152) 



where we have again used z — (t — x/vx)(x/v 2 — t). Note gzigil — 9v29i2* > an d 9§29§2* can 
be obtained from gu.gn by interchanging v\ and v 2 - 

The first thing to observe is the dependence of these quantities on the time difference 
At = t — t'. In the time-ordered case (|151|) the dependence is approximately 1/At, while in 
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the retarded case ( |152|) it is 5 (At). This suggests that in a given configuration both g and g 



R 



are rapidly varying functions of time. The remaining integral in the time-ordered case (|151|) 
can be evaluated for the special case of equal times At = 0, see Eqns. (|C7|) , ( |C8| ) , and flC9|). 
One finds that \g%j(t, x, 0)| 2 is the same as \g^(t, x, 0)| 2 , provided one makes the replacement 
9(t)5(0) i— > l/2ne. Comparing Eqns. ( |146| ) and ( |152| ) we see that the structure of the equal- 
time expressions, \gij(t, x, 0)| 2 and |g^(t, x, 0)| 2 , is similar to the result for D R (t, x, 0), up 
to an infinite prefactor. In particular, the diagonal elements of \g\ 2 and \g R \ 2 have delta- 
functions with exponentially decaying amplitudes, and all elements have a term which decays 
algebraically at large x. 

In comparing the expressions for D R ( |146|) and |g H | 2 ( |1 52|) , one obvious difference is 
the presence of the factor V[l/(tQ — t)] in the density-density correlation function. One 
consequence of this factor is that it makes D R (t, x,0) at fixed x an odd function about 
t — to. We shall now demonstrate that this antisymmetry is present in each realization of 
disorder, not just in the averaged quantity. 

The derivation of this antisymmetry relies on some of the results derived in Appendix [A] 
that allow us to express the Green's function for the disordered 331 sequence in terms of a 
coordinate-ordered exponential. We begin with an expression for the time-ordered, single- 
particle matrix Green's function in the presence of an arbitrary scalar potential £(x), ( |A16| ) 



G^t, x, x') = [ p,e~^ [9(t) n(-uj) - 6(-t) n(uo)} Q~ 1/2 S(x, x'; uo) Q~ 1/2 , (153) 

where n(uj) = [exp(/3uj) + 1] _1 is the usual Fermi distribution function (with /3 the inverse 
temperature), and from Eqns. ( |A5| ) and ( |A10| ): Q = v n H + (\/ix)a x and 



S(x,x'; u) = T y exp J dy 



n -i £(y) 

ojQ a 



(154) 



y/ViV 2 

Eq. (|153|) is analogous to Eq. ( p7|) for the 110 sequence. Both equations express the Green's 
function in terms of a coordinate-ordered exponential. The expression for the 331 sequence 
is more complicated because the term in Eq. (|154|) describing propagation in the absence of 
disorder, uQ -1 , does not commute with the random field term, ^(y)u z . 

Since Q^ 1 = (l/viV2)(v n t. — (\/7r)<r x ), we can factor out an overall phase from S(x, x'\ uj): 



-a H a 



= e iuito s(x,x'; u) 



S(x,x'; u) = e iuJVn{x - x ' )/viV2 T y exp [ -i dy 

\ Jx' \_nv\V2 s/v\V2 jy 

(155) 

Using this result and Eq. ( |1 29|) to transform G$ to g we find 

g(t,x,x') = [Q- 1/2 0} ] [ —e^ t0 ^[6(t)n(-uj) - 6(-t)n(oj)]s(x,x'; u)[Q-^ 2 0]. (156) 
From this equation we find 

g*(t + t, x, x') = -[Q- 1 / 2 0f[ ^e^[6(t + t)n(-u)-9(-t - t)n(u)]s*(x, x 1 ; u) [Q- 1/2 0}\ 

J mi 

(157) 
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From the form of the matrix s (|155|) , one can show: 

s*(x,x'] uj) = a y s(x,x'] u)a y . (158) 
Combining this with the fact that for \t\ < |to| : 

9(t + t)n(-oj) - 6(-t - t)n(u) = 9(t - t)n(-u) - 6(-t + t)n(u>), (159) 
implies that Eq. (|157|) can be rewritten 
g*(t + t,x,x') 

= -[Q- 1/2 0fa y f ^e^[9(t - t)n(-u>) - 6(-t + t)n{u)]s{x, x'; u)a y [Q- l ' 2 0\* 



= -[Q-v<0] o» [[Q- l "0] 1 ) g(t -t,x,x')[Q-^ 2 O] a v [Q-^0] . (160) 
Therefore, if we define the matrix 

C S [Q^O\\' [iQ-MY 1 = ( ° , (161) 

V / \yv 1 /v 2 J 

we have the final result that 

g{t + t,x,x') = -Cg*{t -t,x,x')C T , for \t\ < \t \. (162) 

One can show by an analogous derivation that the corresponding retarded function, 
g R (t,x,x'), obeys exactly the same relation. 

Using the relation between D and g ( 147|) and the fact that det C = 1, Eq. ( |162| ) implies 



for |t| < \t \ 

?R,eD(to + t, x, x') = — 9fte D(t — t, x, x'), QmD(to + t, x, x) = ^smD(t — t, x, x). (163) 

Similarly, using the expression for D R in terms of g and g R ( [L49|) , Eq. ( |162j ) and the corre- 
sponding relation for g R gives 

D R (t + t, x, x) = -D R (t - t, x, x'), (164) 

where we can drop the restriction \t\ < to since for times t outside the interval [0, 2to], D R is 
identically zero from Eq. ( |146| ). We have found that the retarded density-density correlation 
function in a given sample is antisymmetric about the mean arrival time, indicating that 
the point t = t is special, independent of ^ n (^)- This result is true at any temperature. If 
we include disorder in the velocity so that v n is a function of x, the symmetry (|164[) would 
be absent and the principal value singularity in D R would be rounded out, cf. Appendix 0. 

Next we consider the behavior of correlation functions for some simple potentials ^ n (^)- 
We first consider the case of a uniform potential. As we remarked in Section [II B| , the 
correlation functions in the presence of a parallel magnetic field can be reinterpreted as 
those in the presence of a uniform scalar potential difference between the layers £, n (x) = T. 
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Using Eqns. fl40|) and ([TJ) to go from the result in Section [II IB] for the time-ordered single- 
particle Green's function, ([T4|), to the corresponding retarded function, and Eq. (|129[) to 
transform this to the Majorana basis, we find 



g R (t,x) 



2i9(t) 



v n t - -a z ) ^eP x (r,X) + ( -1 - v n a z ) -KeP T (r,X) 

71 I \7T J TV 



VlV 2 

w lV2 TameP{T,X) 



(165) 



The quantity 3?e P(t, X) and its derivatives are evaluated in Appendix O, and from 
Eqns. ( |CT0|) , (|CITD , and (PT2|) we find 



g&(t,x) = -iO(t) 



8{x - fit) - ^B{z)J\ (^j^y/v\v 2 z 



Ix — v 2 t 

V\t — X 



g? 2 (t } x) = -t6(t)—6(z)J 



A/tt 



y/ViV 2 Z 



(166) 



where the J n are standard Bessel functions. Note g^\ = —g\ 2 and g 22 can be obtained from 
gj\ by interchanging v\ and t>2- These functions are plotted in Figs. [II] and |12 . 




FIG. 11. The imaginary part of plotted for v n = 1, x = 10, A/-7T = 0.1, and T = 3. The 
horizontal axis is time measured from the mean arrival time. 
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-10 h 

FIG. 12. The imaginary part of g R 2 plotted for v n = l,x = 10, X/tt = 0.1, and T = 3. The 
horizontal axis is time measured from the mean arrival time. 



We see that the delta-functions present in the diagonal elements of g R in the case of 
zero potential (Eq. ( |150| ) evaluated at = 0) remain in the case of a uniform potential. 
In addition to these delta-functions we find oscillatory terms in all elements of g R (t, x) 
considered as functions of t at fixed x. Since the quantity z — (t — x/vi)(x/v 2 — t) is 
maximal at t = to, we see that the frequency of the oscillations in g R is minimal at the mean 
arrival time and increases as we approach the extremal arrival times. This feature is evident 



in Figs. O and p. 



Finally, we consider the case of a scalar potential made up of isolated impurities located 
at the points {y m } with strengths {q m }, i.e., we take the potential to be 

in{x) = ^q m 5(x - y m ). (167) 

m 

Note that the white-noise potential used in the previous calculations can be approached by 
the form given in Eq. (|167|) if we take the number of impurities to infinity and the q m to 
be random variables. Using Eq. (|126|) , one can compute g for this potential by an exact 
summation of the perturbation expansion. Once again, it is the chirality of g^ that makes 
the calculation tractable. 

If there are N impurities between x and x' , then the nonvanishing terms in the per- 
turbation expansion of gij are in one-to-one correspondence with the set of all iV-tuples of 
nonnegative integers. For example, if x' < y± < y<i < . . . < i/n < x, then for t > the 
nonvanishing terms correspond to propagation from x' to y\ followed by scattering n\ times 
off impurity q%, followed by propagation from y\ to y<i and scattering 712 times off impurity 
q 2 , etc., where n m are nonnegative integers. If n m is even then the velocity of the particle is 
unchanged by the scattering [i.e., the internal lines on either side of the impurity are either 
both gfi or both g^), while if n m is odd then the velocity of the particle is changed by the 
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scattering (i.e., the internal lines on either side of the impurity are g$ and #22 )■ Therefore, 
from the parity of each element n m of the iV-tuple we know which internal propagators are 
g$ and which are #22 , an d we can define a corresponding arrival time T. This arrival time is 
given by T = Xi/v\+ X 2 /v 2l where Xj is the total distance the particle travels with velocity 
i>j. From these considerations we can conclude that the general form of the time-ordered 
and retarded Green's function for the case of isolated impurities is: 



gij(t,x,x') =J2 L 



1 



g%(t,x,x') = -ie(t)J2^S(Ti j -t), (168) 

k 

where the arrival times depend on the end points x, x' and the location of all impurities 
y m located between the end points, and the coefficients L % 1 depend on the impurity strengths 

Qm ■ 

In calculating g^ with N impurities between x' and x, these impurities divide the interval 
[x 1 , x] into N+l segments. Since we are computing gij, the first segment must be associated 
with a factor of g^ and the last segment with a factor of g^ . Since on the N — 1 remaining 

segments we can have either g^± or g^ depending on the parity of the n m , there are in 
general 2 N ~ l arrival times Tj? , and thus the number of terms in the sum in Eq. ( |168| ) grows 
exponentially with the number of intervening impurities. 

For a given arrival time T^- 7 , the parities of the Tl m S £1X6 determined and the corresponding 
coefficient L % 1 in Eq. ( |168| ) is given by 



JV 

L k = C ij II r m, Pm , (169) 
m=l 

where Cu = 1/fj, c i2 = ij y / viv 2 = — c 2 \ are overall coefficients, p m = ra m mod2 is the parity 
of n m , and 



r m,0 = 2^ ~1==^F 



2n -| 



oVv 7 ^/ l + ?m/ui«2 

f ( i( lm V n+1 iq m /VviV2 n „„, 
r m>1 = 2^ — = = v , 170 

where we have assumed \q% n /viV2\ < 1 for all m. The coefficient r m>0 (i) is the sum of the 
amplitudes for scattering off the mth impurity an even (odd) number of times. 

For a given realization of the potential in Eq. ( |16?1 ), the Green's functions in Eq. ( |168| ). 



considered as functions of time, have an exponentially large number of singularities. How- 
ever, the disorder-averaged quantities, g (|128|) and g R (|150[) , are very simple functions of 



time. This is readily understood if we approximate the white-noise potential by taking the 
impurity strengths q m in Eq. ( |167| ) to be independent, identically-distributed random vari- 
ables with zero means. Then from Eq. ( |1 70| ) we see r^J = for all m. Thus the only 
nonvanishing L % 1 involves an even number of scatterings off each impurity, and hence there 
is only one singularity, corresponding to the particle never changing its velocity: 
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9ij(t, x, 0) = dij— — — (r mfi ) N . (171) 

2tt(x — x' — vd + ie t ) 



From Eq. ( 170 ) we see that r^o < 1, and since N, the number of impurities between x' and 
x, is proportional to \x — x'\, the last factor in Eq. (|171|) reproduces the exponential decay 
present in Eq. ( |128| ). 

We can understand several things from the general form of the Green's functions given in 
Eq. ( |168j ). First note that the fact that g R for the case of isolated delta-function impurities 



is a sum of delta-functions in time ( 168 ) is consistent with the S(t — t') factor present in 



gfj(t,x,0)gfj*(t',x,0) ( |152j ). One can similarly show that the form of g in Eq. ( |168| ) is 
consistent with the — t') behavior of g%j{t, x, 0)g*j(t', x, 0) ( |151| ). 

Next we suppose that the N impurities between x' and x are evenly spaced. In the limit 



N — > oo it should not make a difference. In this case, Eq. ( |168| ) still holds but the quantities 
VI and T]? are determined by a different set of rules, so we shall distinguish them by a 
tilde. The number of arrival times T^ 3 is reduced from 2 N ~ 1 down to N. This is because 
the interval \x' — x\ is now divided into N + 1 segments of equal length. The velocity on 
the first and last segments is fixed by the indices on g^ and the N arrival times can then 
be specified by the number k of the N — 1 remaining segments which have velocity v 2 . For 
example, for gn the arrival times are 

~n / k N + 1 - k\ \x-x'\ . 

T "=U + ^— y-WTT- (172) 

where k = 0, 1, . . . (N — 1). The amplitude L % 1 is the sum of ( A? ^ 1 ) of the L^'s given in 

Eq. ( |169| ) . In the limit of large N the number of terms which contribute to L% is a Gaussian 
distribution in k peaked at k = (N — l)/2. The arrival time corresponding to this amplitude 



is 



~ n /JV-1 N + 3\\x-x> 

-if/V-lV? — ~ ^ r 



~ \ 2v 2 ^ 2v x ) N+l 
v n \x-x'\ = ^ 

V\V 2 

Thus we find that the mean arrival time to emerges in this model potential because the 
number of terms that contribute to the singularity at is maximal for = to. 

All of the considerations in this section lead us to the conclusion that D R (t,x,x') in 
a typical configuration is given by D R (t,x,x'), Eq. ( |146j ) , with the factor of V[l/(t Q — t)] 
replaced by a function f(t,x,x'), which is a rapidly fluctuating function of time, antisym- 
metric about the point t = to, and whose amplitude grows as t approaches to. The claim 
that the general structure of D R in a typical configuration is captured by D R is supported 
by the relation between D R and g, g R , Eq. ( |149|) , and the fact that the second moments \g\ 2 
and \g R \ 2 , Eq. (|152 ), have the same structure as D R , without the principal value factor. We 



expect that /(£, rapidly fluctuating function of time based on the dependence of 



g ij (t,x,0)g* j (t',x,0) and g R (t,x,0)g R ^(t',x,0) on (£-£'), E q ns - (EUD and (HI), and the 



behavior of g and g R for the case of isolated delta-function impurities, Eq. (|168|) . However 
complicated f(t,x,x') is, we know it must be antisymmetric about t = t by Eq. ( |164| ). The 
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claim that the amplitude of f(t,x,x') approaches a maximum at t = to is supported by 
several results. First, recall that g R in a constant potential T, Eq. (|166|) , is oscillatory with 
a frequency which is proportional to T and which is minimal at the mean arrival time. This 
suggests that in a potential which varies with x there will be less cancelation near the mean 
arrival time than near the extremal arrival times. Second, there is the observation that 
in the model of equally-spaced isolated impurities the number of terms which contribute to 
each singularity is maximal for the singularity at the mean arrival time. Finally, we have the 
fact that averaging produces a function singular at t — t Q . These conclusions are supported 
by our numerical simulations, see Section |TJ. 



D. Pfaffian Edge 

In this section we consider the edge theory of the Pfaffian state, concentrating on the 
form the retarded density response function in the presence of a finite temperature and 
disorder. 

The edge theory of the Pfaffian state contains a c = 1/2 minimal model CFT in addition 



to the usual (c = 1) chiral boson |TJ|. Recall that the c = 1/2 minimal model has three 
primary fields, which we denote as J, x, and a, and whose scaling dimensions are 0, 1/2, and 
1/16, respectively. The x fi e ld is identical to a single chiral Majorana fermion. The chiral 
boson in the edge theory, (p, can be shown to have a compactification radius of R = l/v^8, 
the same radius as the charged mode of the 331 state, see I. In addition to the c = 1/2 
primary fields and their descendents, the operator content of the edge theory includes the 
primary field d x ip and the vertex operators e lkip ^, where k £ Z. The electron creation 



operator is 24 



:&{x):=—x(x)e+ /s «'\ (174) 
lira 

which has a scaling dimension of 3/2. 

We now determine the electric charge density operator, p(x), for the Pfaffian edge. One 
requirement for this operator is that the electron operator has a unit charge with respect to 
it 

[p(x), ¥(x')\ = 8(x - x')^(x'). (175) 
One candidate which satisfies this relation is 

P(x) = ^pf*^ x )- ( 176 ) 

The operator on the r.h.s. of this equation has dimension one. The only other dimension 
one operators present in the edge theory are e ±iv/ ^. However, if we were to add to our 
definition of p(x) some non-zero multiple of the Hermitian combination (e 4 "^* 3 + e~ iv/ ^), 
we would find that the condition ( |175[) is violated. Any other higher- dimension operators 
which could be added to p{x) would necessarily be multiplied by explicit powers of the short 
distance cutoff a, and will therefore vanish in the continuum limit. The Pfaffian state has 
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no interlayer dynamics. Each electron is in a state symmetric between the layers, under the 
assumption that the symmetric-antisymmetric splitting is large. Hence, there is an energy 
gap for any process which excites the layers independently. 

Having determined the charge density operator for the Pfaffian edge theory, we can imme- 
diately write down the retarded density two-point correlation function at zero temperature 
for the clean system 



■—6(t)6'(x-v v t). 

47T 



(177) 



From the results in Section [111 (J| we know that the only RG non-irrelevant disorder is 
scalar potential disorder. This does not modify V§ f at all, just as in the case of the charged 
mode for both the 110 and 331 sequences. From the discussion in Section |11 D| we know that 
at a finite temperature this correlation function is also unchanged. We find that Vp { (t,x) 
has a signal at a single velocity even at a finite temperature and in the presence of a non-zero 
scalar potential. 



E. Numerics 

Here we briefly discuss the numerics for the neutral mode contribution to the retarded 
density response function of the 331 edge. The matrices S(x, 0; uj) were computed using a 
discretized version of Eq. ( |154j ). Using Eq. ( |153| ) the single-particle Green's function was 
then found by a fast Fourier transform (FFT) algorithm. Finally, we used the relation 
between the single-particle Green's functions and the density response function (148[) to 



calculate D R , and integrated over time to find the neutral mode contribution to the signal 
in Eq. (M). 
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APPENDIX A: DISORDERED 331 SEQUENCE VIA THE SPIN- ANALOGY 

In this appendix we present an alternative method for obtaining exact disorder-averaged 
correlation functions for the 331 sequence. It is similar to the method used in our discussion 
of the 110 sequence in Section [HI C| , with some technical complications. The procedure is 
based on the fact that we can explicitly construct solutions to the Heisenberg equations of 
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motion for each realization of disorder. We shall ignore the charged mode throughout the 
discussion. 

The 331 Hamiltonian including a disordered scalar potential ( |123| ) expressed in terms of 
the Dirac field is 



L/2 
-L/2 



dx 



A 



2tt 



(Al) 



where in this appendix we suppress the subscript on £ and take periodic boundary conditions 
ip(x + L) = ip(x). The Heisenberg equations of motion for the field operator if)(t,x) and its 
Hermitian conjugate are: 



[d t + v n d x + if (a;)] ip(t, x) + {X/ir)d x ^(t, x) 
[d t + v Q d x - if (a;)] ip^t, x) + (X/Tr)d x i/j(t, x) 



0. 
0. 



(A2) 



The anomalous (i.e., fermion- number non-conserving) terms in Tip couple the equations of 
motion for if) and ^ and therefore we must expand the field in terms of both creation and 
annihilation operators 



x e 



-iuJnt 



a n + B n (x)e iuJnt al 



(A3) 



where by assumption a) n are canonical Fermi operators which create exact single-particle 
eigenstates ofHo with energies u> n . Substituting this expansion into the equations of motion 
leads to the following matrix differential equation 



iQd x f n {x) = [uj n t - £(x)a z ]f n (x), 



where 



Q 



( V n A/tt 

VA/tt v n 



A n (x) 
B*(x) 



(A4) 



(A5) 



Note that the Hamiltonian (|A1|) has a particle-hole symmetry under which ijj <-> 
and £ i — ► — £. In terms of the two-component wavefunction f n {x), this implies that if f n is a 
solution to Eq. (|A4j) with energy u n then f n = a x f* is a solution with energy — u n . Assuming 
all uj n 7^ 0, we can enumerate the functions f n (x) in such a way that uj_ n = —oj n and u n > 
for n > 0. This implies /_„ = /„, from which we find A n = B_ n , an indication that some 
double counting may be present. Indeed, the particle-hole symmetry only interchanges the 
two equations (|A~2] ); it should not generate new solutions. This double counting can be 
removed if we define a n + aL n t— > a n , and write, instead of Eq. (|A3|), 



(ip{t,x) 



x e 



-ILOnt 



To obtain the solutions to Eq. (SI) we define the rescaled wavefunctions 

Cn Q ^ fm 



(A6) 



(A7) 
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in terms of which the differential equation ( |A4|) becomes 



d x ( n (x) 



u n Q 1 - 



fa) t 



We can write the solutions to this equation as 



with a coordinate-ordered exponential 

S(x, x'; u) = T y exp \i I dy 



uoQ 



-i 



*(y) 

( 



(A* 



(A9) 



(A10) 



for x > x 1 ', and the Hermitian conjugate S(x,x'; to) = S\x' ,x; u) for x < x' . The boundary 
conditions on the Fermi field imply ( n (x + L) = ( n {%), which in turn means that the allowed 
energies u n are determined by finding those energies for which the matrix S(L, 0; u) has a 
unit eigenvalue with the corresponding eigenvector taken to be Cn(0), 



[S{L,0; w) - 1]C„(0) = 0. 



(All) 



The orthogonality of solutions f n (x) for different values of u> n is guaranteed by the fact 
that the differential equation ( |A4[ ) is self-conjugate, while their normalization must be de- 
manded explicitly, 



L/2 
-L/2 



dxfl(x)f n (x) = 1. 



This can be rewritten with the help of Eqns. ( |A7| ) and ( JAT5| ) as 



VlV 2 



L^C t (0)AT(^,OC(0) 



(A12) 



(A13) 



A 



7lV n L 



L/2 
-L/2 



dxS j {x,0; u)a x S{x,0; 



to 



Usually it is the wavefunction normalization that makes the disorder calculations so difficult. 
Note, however, that the integration in the second term of the normalization matrix J\f(u, £) 
is extended over the entire length of the sample, which makes it a self-averaging object. As 
in Section [111 Q , the matrix S(x, x'\ uS) ( |A10| ) can be interpreted as the evolution operator in 
a fictitious time y, for a spin precessing under the influence of a constant magnetic field in 
the x-direction (due to the off-diagonal terms of the matrix Q~ l ), and a random field ~ ^(x) 
along the z-direction. The integral in Eq. ( |A13| ) is the x-component of the spin averaged over 
a "time" L. In the presence of a non-vanishing disorder, the initial orientation is forgotten 
after a finite distance, and the second term in the normalization matrix ( |A13| ) disappears 
in the thermodynamic limit, L — ► oo. This happens with probability one for any realization 
of disorder. Physically, this simplification is related to the fact that in a chiral system 
localization does not happen; each particle explores the entire circumference of the sample. 
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Let us now consider the single-particle Green's function of the neutral mode fermion, Q$ 
[Eq. (fflD l. Using Eq. (|A6|), this can be written 



G4t,x,x') = -iY,fn{x)ft{x')e- i ^ t [B{t)n{-u n )-e{-t) 



nlu n ) 



(A14) 



where n(oS) = [exp((3u) + 1] 1 is the usual Fermi distribution function. Using Eqns. 
and (|A13|) , with Af(u,£) = 11, we obtain 



Gi,(t,x,a/) = -i ( -^J Y.< 



-iu) n t 



xQ- 1/2 



[0(t)n(-u n )-6(-t)n{u n )} 



S(x, 0, u n ) C " (0)C " (0) gt(^ o, u n ) 



Q 



-1/2 



(A15) 



where Cn(0) obey the eigenvalue equation (|A11[ ). To take the thermodynamic limit, we need 
to set the system size to infinity, keeping other parameters (temperature, disorder, distance 
\x — x'\, etc.) finite. Effectively, this implies that we can select an energy interval AE, 
"infinitesimal" on a scale defined by these finite quantities, and yet containing a macroscopic 
number of energy levels, such that the averaging over the states within this interval gives 



CnCn) rAF 



1 



The value of the average and the existence of such an interval follows from the fact that 
for any uj\ ^ u 2 , the spin-rotation matrices S(L, 0; UJ\^) become entirely uncorrelated for 
a sufficiently large L, or, equivalently, the relative rotation matrix S^(L, 0; oj\) S(L,0; u> 2 ) 
entirely forgets the initial direction. 

Performing the averaging over the eigenstates within such an interval, we obtain for the 
correlation function ([A15Q , 



G$(t,x,x') = -i 



dui 
— < 
2vr 



-iujt 



[9(t) n(-u) - 6{-t) n(uo)} Q~ 1/2 S(x, x'\ u) Q~ 1/2 , (A16) 



where the summation was replaced by an integration using the "clean" single-particle total 
density of states p = (ff 1 +v 2 ~ 1 )L = 2Lt> n /t>it> 2 , which cannot be modified by disorder. Note 
that Eq. (|A16 ) does not contain a disorder average; it is an expression valid for any given 
realization of disorder (or even in the limit of no disorder, as long as this limit is taken after 
the thermodynamic limit). For example, we checked that Eq. ( |A16| ) with £(x) = const, 
reproduces Eq. (0), which was derived by more conventional methods. 



From the definition of S(x,x'; u), Eq. ( |A10| ), and the disorder averaging procedure used 
previously, see Eq. (|113|) , we find 



e -A i \x-x'\/2v 1 V2 e Mx-x')Q- 1 



S(x, x'\ io) 

With the help of Eqns. (|70|) and ( |A16|) this gives, in the zero temperature limit, 



(A17) 



G(t,x,0) 



[27r(x - v c t + ie t )} 



1 

m— 1 2 



+ 



1 - a a 



2ti(x — Vit + ie t ) 2ir(x — v 2 t + ie t ) 



-A^\x\/2viV2 



(A18) 
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This is in exact agreement with Eq. ( |130| ) of Section [III C[ We have checked that the other 
disorder-averaged correlation functions for the 331 sequence discussed in Section [III C| can 
also be reproduced using the method described in this appendix. 



APPENDIX B: RANDOM TUNNELING FOR THE 331 SEQUENCE 

In this appendix we illustrate the effect of an RG irrelevant random perturbation by 
analyzing the neutral mode of the 331 bilayer in the presence of velocity and tunneling 
disorder. Specifically, we assume that both the neutral mode velocity v n (x) and the tunneling 
amplitude X(x) in Eq. ( [All ) are coordinate-dependent, in such a fashion that the system 
remains chiral, vi^x) = v n (x) ±X(x)/tt > for all x. The introduction of such a coordinate 
dependence requires only a slight modification of the Hamiltonian (|A1|) . Specifically, the 
first term in the Hamiltonian density must be replaced as 

-iv n ^d x ij) -> -~ \y n (x) il) ] d x ip + ip ] d x [f n (x)^]} . 

The arguments in Appendix ^ can then be repeated with little modification and we obtain, 
in place of Eqns. ( |AT5| ), (|&I0D , 

G^t, x, x') = [ ^e~ iU [6(t) n(-uj) - 6{-t) n{uj)} Q- 1/2 (x) S(x, x'; lu) Q- 1/2 (x'), (Bl) 

i ( f x r iA £{y) 

S(x 1 a/;u) = T y exp[i dy uQ~ (y) - u{y)cr z , u(y) = 1/2 1/2 (B2) 

Again, this expression is valid for any given configuration of disordered v n (x), X(x), and 

The requirement vi^ix) > is equivalent to essentially non-Gaussian disorder, and the 
disorder averaging is generally non-trivial. This, however, is greatly simplified in the absence 
of potential disorder, = 0. In this case the remaining matrices Q~ x {y) in the exponential 
commute with one another for all y, the coordinate ordering (T y ) can be omitted, and the 
disorder averaging can be performed directly. The structure of the expression is most evident 
after the unitary transformation (|1 29|) to the Majorana fermion representation, 



gij(t, x, x') = O^G^{t, x, x') O 

-J^n 6 m n( ~"-° " ^ <UJ)] v}'\x)v] l2 {x>) ' { ) 

where Ti(x, x') = f®,dy/vi(y) is the time it takes for the iih mode to travel from x' to x. 
Clearly, the structure of the correlation function in a given configuration of disorder does 
not change; at zero temperature we obtain [cf. Eqns. (|(J), ( |13U| )] 



G^(t,x,x') = — 



H + a x II - a 1 

+ 



_Vx 2 {x) v I (x 1 ) (ti(x,x') — t) v^ 2 (x) v^ 2 (x') (t 2 (x,x') — t) 



The disorder averaging can be performed for weak disorder if we notice that the velocity 
fluctuations along the entire path contribute to the arrival times Ti(x,x')] these quantities 
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acquire nearly Gaussian distributions at sufficiently large distances (compared to the disorder 
correlation length, l c -C \x — x'\). If we ignore small multiplicative corrections near the ends 
of the interval, we then find 



G^(t,x } 0) = — 



-1/2 



a + 0- 



Di 



DAx) = tHx) 



where Tj = Tj(x) = Ti(x) — t is the time elapsed from the arrival of the ith peak, D ri = 
T 2 (x) is the corresponding dispersion, and 

F (T)= r^e^-" 2 ' 2 , (B4) 



h 2iri 

where we have assumed t > 0. As illustrated in Fig. [13], at large values of the argument, 
this function approaches asymptotically the clean single-particle Green's function, F(T) = 
(2ttT)^ 1 , \T\ ^> 1. Although the perturbation is RG irrelevant, the form of the Green's 
function is modified in the vicinity of the singularities. 
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FIG. 13. The real and imaginary parts of the universal function F(T) [Eq. ( p34j) 1 , which de- 
scribes the shape of the peaks of the averaged Green's function for the 331 double layer. Dotted 
line shows the real part of the Green's function in the absence of disorder. 



For a weak (w v n ) Gaussian disorder with a finite correlation length, 
X(x) = X + 5X(x), (5X(x)5X(y)) = w 2 f(x-y) /(0) = 1, 
assuming w 2 C we obtain, to leading order in the weak disorder expansion, 



n {x) = x/v [1 + w 2 /7r 2 vf + 0(w A /vt)], 
where Di(x) = w 2 xl c /vf + . . ., and the disorder correlation length, 



dx f(x) 



was assumed to be short compared with the overall distance, l c <^ \x\ 
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APPENDIX C: EVALUATION OF INTEGRALS 



In this appendix we evaluate the integrals needed in the main text. The basic integral 
is of the form 



oo 



A , Y s i ^ T sinh(A(a;)X) 
j-oo 2ir A-yu) 



where A (a;) = \/ a 2 — b 2 uo 2 , and r,X,a, and b are real parameters. At large frequencies the 
integrand is of the form 1/uo times an oscillatory function, and therefore converges without 
a regulator. Also, since the expansion of the integrand in powers of A(u;) contains only 
even powers, no branch cut is needed and the integrand is therefore an analytic function for 
finite lo. For large lo the integrand contains the factors e lul( - T±bX \ Thus for \bX\ < |r|, the 
integrand is exponentially small on one side of the real axis and we can therefore close the 
contour on that side and find A = 0, hence 

A(t,X) ot 6[(bX) 2 - t 2 }. (C2) 

To evaluate the integral in the region for which it is nonzero, we break up the integral 
into two terms: 

r duo p a H x r du f -H^)x 
v ' ; Jci 2tt 2A(w) Jc 2 2tt 2A(uj) ' v ; 

where now we must introduce a branch cut, which we take to run along the real to axis from 
—a/b to a/b. The contours C\ and C2 are shown in Fig. [14]. Assuming X > 0, the contours 
have been chosen so that the integrals in Eq. ( |C3| ) are separately convergent. The integrand 
in the C2 integral is exponentially small in the upper half plane and thus can be closed there 
to give zero. The integrand in the C\ integral is exponentially small in the lower half plane 
and can be closed there and contracted to run around the branch cut. We use the change 
of variables u = {a/b) simp in the integral around the cut to arrive at 

A(r,X) = e\{bXf - T 2 ]\ 1^ ^ e a[i(T/6)«m V +Xc«d (q 4 ) 

2b Jo 2tt 




FIG. 14. The complex u plane showing the branch cut along the real lo axis and the contours 
C\ and C2. 
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We next note that we can write 



. T „ /„„ /T N 2 



i— simp + X cos (f = d X 2 — y—j cos(<£> — i(po), (C5) 



where the real parameter ip$ is defined by cosh<y2o = X/JX 2 — (r/b) 2 . After performing 



this substitution in Eq. (CJ ) and noting that the integrand is a periodic function in <p with 



period 2ir } we can perform a final change of variables 6 = <p — ic^o to find 
A(r,X) = 9[{bXf - r 2 ]l ^v^M^cos, = ^ [(feX)2 _ ^ 



2 _ r 2 



(C6) 



where Jo is a Bessel function of imaginary argument. 

Using a = Af, b = \/it,t = v n X — i and differentiating Eq. (|C6| ) with respect to r and 
X gives the following results 

:t^ sj ^-^H^A (c7) 

00 ^ e Ma-f) cosh(A(^)X) = -sgn(t) [5(x/v 1 -t)+ 6(x/v 2 - t)] + -^-6(z) 
-00 2w 2 2viV2 

x ^ (^) • < C8 > 

00 do; ,• ,„ Y sinh(A(co>)X) 1 , . r _. . , . 

e w(tJ " f) iw ' v ; ; = -— — sgn(t) pfa/ui - i) - 5(x/v 2 - i)] 



2tt " A(w) 2A/vr 



A? r W" nX/ T 2 ~M #^l, (eg) 



2(A/tt) 2 v 7 ^ 1 VA/tt 

where 2 = (i — x/vi){x/v 2 — t). 

Finally, if we analytically continue the above results to the case of A^ = i^/v^U^T, where 



T is real, and write k(lu) = Jv\v 2 T 2 + (X/tt) 2 uj 2 , we find 



-00 2tt k(u) 2(A/tt) v ' \X/n 

dui 
2^ 



z)Jo h-Vw , (CIO) 



v <Lj Mv n x-t) cos ( K ( tu ) X ) = Lgn(t) [5(x/v 1 -t) + 5(x/v 2 - t)] - ^-9{z) 



x X J x \^—^/ViV2z \ , (Cll) 

-e"^^ V ; ( ) = -— — sgn (t) [S (x vi -t)- S(x v 2 - t)} 
-00 2ir k{u>) 2A/7T 

T8(z) (VnX-V^t) T ( T y /v 1 v 2 z \ 
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